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ABSTRACT 

Experiments  were  performed  in  which  the  velocity  pro¬ 
files  and  decay  and  spreading  properties  of  free  underexpanded 
Jets  of  cold  air  were  measured.  Stagnation  point  heat  trans¬ 
fer  parameters  for  impingement  of  these  jets  on  various 
surface  shapes  were  evaluated  and  correlations  made  with  the 
free  Jet  data.  Pressure  distribution  and  photographic 
studies  of  the  free  and  impinging  Jets  revealed  that  an 
unusual  separated  flow  phenomenon  can  exist  under  certain 
impingement  conditions.  Problems  associated  with  changes  in 
Jet  stability  and  the  effects  of  interference  due  to  geometric 
arrangement  of  the  apparatus  were  considered  in  a  qualitative 
way. 
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1.  INTRODUCTION 


The  problem  of  estimating  the  stagnation  point  heat 
transfer  from  a  jet  of  hot  gas  to  a  cold  surface  upon  which 
it  impinges  can  be  treated  by  applying  the  usual  techniques 
of  stagnation  point  heat  transfer  computation.  These  solu¬ 
tions  are  generally  expressed  with  the  stagnation  point 
velocity  gradient  (due/dr)r_Q  as  a  parameter  (see,  e.g. 

[l]).  For  the  case  of  the  impinging  Jet,  this  parameter 
can  be  readily  determined  from  the  pressure  distribution 
measured  on  the  impingement  surface .  In  order  to  provide 
suitable  data  for  such  determinations,  a  broad  experimental 
study  was  made  of  a  number  of  different  Jet  impingement 
cases.  This  study  included  the  impingement  of  an  unheated, 
turbulent,  axially  symmetric  Jet  of  air  on  several  different 
surface  shapes.  The  Jet,  which  issued  from  a  convergent 
nozzle,  was  run  at  several  pressure  ratios,  both  subsonic 
and  sonic  (underexpanded),  and  the  impingement  distance  was 
varied  through  a  wide  range.  The  surface  shapes  used  were 
a  convex  hemisphere,  a  flat  plate,  a  concave  hemisphere,  and 
a  shallow  cylindrical  cup. 

In  addition  to  the  stagnation  point  radial  velocity 
gradient  itself,  the  correlation  of  this  gradient  with 
measured  free  Jet  properties  was  evaluated.  One  correla¬ 
tion  was  based  on  properties  of  the  free  jet  at  the  nozzle 
exit,  and  another  was  based  on  local  properties  of  the  free 
Jet  at  the  same  location  as  the  impingement  surface.  These 
parameters  are  discussed  in  3.2.3. 

Although  the  primary  characteristics  of  turbulent  free 
jets  (axial  decay,  radial  spread,  etc.)  are  well  known, 
several  interrogated  secondary  factors  which  influence  their 
detailed  structure  should  be  considered  if  one  is  to  corre¬ 
late  data  resulting  from  tests  made  with  different  nozzles 
and  supporting  structures.  In  the  present  context,  in  which 
It  is  assumed  that  solid  boundaries  may  exist  within  a  few 
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nozzle  diameters  of  the  Jet  source,  the  most  important  of 
these  factors  are  thought  to  be  (a)  a  dynamic  instability 
effect  and  (b)  a  blockage  or  interference  effect.  The 
dynamic  instability  referred  to  is  characterized  by  a 
primarily  lateral  oscillation  or  "flapping"  of  the  entire 
Jet  flow  field,  and  is  to  be  distinguished  from  the  shear 
induced  instabilities  of  the  turbulent  mixing  flow  field 
which  are,  in  general,  of  higher  frequency  and  lower  ampli¬ 
tude.  The  blockage  effect  is  distinguished  by  changes  in 
core  length  and  spreading  rate  due  to  the  presence  of 
obstructions  in  the  entrained  flow  field.  Since  the  exter¬ 
nal  configuration  of  the  nozzle  itself  can  thus  be  a  factor, 
it  is  important  to  consider  this  effect  in  determining  Jet 
flow  field  momentum  balances.  A  Jet  with  a  coaxial  free 
stream  velocity  superimposed  about  it  is,  in  a  sense,  equiva¬ 
lent  to  a  free  Jet  with  a  prescribed  entrainment  flow. 

Studies  which  reveal  the  changes  in  core  length  and  spread¬ 
ing  for  different  ratios  of  Jet  to  fi*ee  stream  velocity  are 
thus  indicative  of  some  of  the  effects  to  be  expected  due  to 
blockage.  In  addition  to  the  two  secondary  effects  Just 
described,  there  is  an  interdependence  of  the  nature  of  the 
instability  upon  structural  interference  through,  for  example, 
the  reflection  and/or  excitation  of  acoustic  disturbances. 

The  aforementioned  effects  and  their  qualitative  influence 
on  the  measurements  are  discussed  more  fully  in  2.2.4. 

In  recognition  of  the  importance  of  the  above  factors, 
a  thorough  experimental  survey  was  made  of  the  same  free 
Jets  that  were  used  in  the  impingement  studies  using  the 
identical  test  setup.  In  this  way,  it  was  hoped  that  mere 
meaningful  correlations  of  impingement  heat  transfer  param¬ 
eters  could  be  obtained. 

Also  of  interest,  in  the  over-all  program,  was  the 
evaluation  of  effects  due  to  the  shock  structure  present  in 
tue  core  of  underexpanded  Jets.  Decay  and  spreading  parameters 
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were  determined  for  a  number  of  combinations  of  jet  pressure 
ratio  (subsonic  and  under expanded)  and  axial  location.  The 
results  of  this  free  Jet  study  as  well  as  a  general  discus¬ 
sion  of  free  jet  structure  are  presented  in  Section  2.  The 
impingement  studies  and  the  correlations  based  on  the  free 
Jet  data  are  treated  in  Section  3. 
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2.  FREE  JET  STUDIES 

trueture  of  the  turbulent  axially  symmetric  free  .let. 

Although  an  extensive  literature  on  the  structure  of 
turbulent  free  Jets  exists  (see,  e.g.  [2]),  relatively 
little  quantitative  information  is  available  on  decay  and 
spreading  behavior  for  underexpanded  cases  [3, ^,5, 17, 23] , 
Analytical  and  semi -empirical  methods  for  determining  these 
characteristics  have  usually  been  restricted  to  either  sub¬ 
sonic  or  properly  expanded  supersonic  Jets.  In  the  present 
study,  however,  it  is  of  interest  to  determine  not  only  the 
local  details  of  the  flow  within  the  jet  core,  but  the 
effects  of  underexpansion  on  decay  and  spreading  rates  in 
the  downstream  portions  of  the  Jet.  Only  recently  has 
interest  in  the  characteristics  of  rocket  exhaust  plumes  at 
high  altitudes  and  in  space  spurred  efforts  to  understand 
such  structural  details.  Some  specific  problems  which  have 
been  studied  involve  pressure,  thermal,  and  shock  inter¬ 
ference  effects  on  adjacent  structure,  as  well  as  vehicle 
stability  effects,  and  the  blocking  of  communication  signals 
due  to  ionization  radiation  in  the  plume.  The  emphasis, 
however,  has  been  on  the  determination  of  initial  spreading 
boundaries  of  the  plume  and  the  strength  and  location  of 
the  initial  shock  structure  [6-17,21,22,24,25,26]  rather 
than  decay  processes.  Other  recent  studies  have  been  con¬ 
cerned  with  jet  applications  in  the  fluid  amplifier  field 
[l8]  and  the  interrelation  between  underexpanded  Jet  stabili¬ 
ty  and  associated  sound  generation  phenomena  [l9l. 

The  general  structural  features  of  turbulent  free  jet 
flows  are  well  known.  If  we  consider  the  flow  issuing  from 
a  simple,  circular,  convergent  nozzle*,  three  major  variations 

^Because  of  the  specific  interest  in  underexpanded  jets, 
the  convergent  nozzle  was  considered  to  be  well  suited  to  a 
study  of  basic  effects  since  there  is  no  dependence  of  the 
degree  of  underexpansion  on  area  ratio  or  nozzle  divergence 
angle.  The  following  discussion  is  thus  limited  to  convergent 
nozzle  flows,  and  the  description  of  certain  aspects  such  as 
shock  formation  is  not  to  be  considered  general. 
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of  the  flow  pattern  are  possible,  depending  upon  the  pres¬ 
sure  ratio  through  the  nozzle.  (Although  a  properly  expanded 
Mach  1  Jet  can  exist  in  principle,  it  is  not  treated  in  the 
present  discussion.)  The  idealized  structural  features  of 
each  of  these  variations  as  well  as  the  nomenclature  and 
symbols  used  to  describe  these  Jets  throughout  the  remainder 
of  this  report  are  shown  in  Figure  1.  The  pressure  ratio 
values  given  are  for  air.  A  typical  sc.hlieren  picture  of 
each  Jet  type  is  shown  in  Figure  2. 

2.1.1.  Subsonic  jet.  A  region  of  turbulent  mixing  between 
jet  and  ambient  fluid  begins  to  form  a  short  distance  aw&y 
from  the  nozzle  lip.  Radial  diffusion  or  spreading  of  this 
region  continues  both  inward  and  outward  as  the  distance 
downstream  increases  until  finally  the  inward  diffusion 
reaches  the  Jet  axis.  At  this  point,  the  "potential"  core 
ends,  but  the  outward  diffusion  and  entrainment  of  ambient 
air  continue.  After  a  so-called  "transition"  region,  at  a 
point  somewhat  farther  downstream,  the  decay  of  axial  veloc¬ 
ity  on  the  center  line  and  the  radial  spread  of  the  velocity 
profile  behave  in  a  manner  consistant  with  the  self -similarity 
of  velocity  profiles  from  that  point  on.  The  Jet  is  now  said 
to  be  fully  developed.  For  air,  the  Jet  will  be  subsonic  for 
isentropic  pressure  ratios  1  >  p^/p^  >  .528.  Throughout 
this  range  it  can  be  assumed  that  p,/p  =  1. 

It  should  be  noted  that  an  ideally  expanded  supersonic 
Jet,  which  contains  no  shock  waves,  has  essentially  the  same 
structure  as  the  subsonic  jet.  Effects  due  to  compressibility, 
however,  become  much  more  important  in  determining  the  core 
length  and  decay. 


^"ideally  expanded"  refers  to  a  properly  expanded  Jet 
issuing  from  a  nozzle  with  zero  exit  divergence  angle,  i.e. 
with  parallel  flow  at  the  exit. 


Sta .  1 


1,1  <  Pl/Poo  <  2 


'////.  Intercepting  ^Reflected  oblique  shock 
■^4///  oblique  shock 


Highly  underexpanded  jet 


°26  >  Poc/Psc  >  0 

2  <  Pi/Poo  <  00 

Figure  1,  Three  major  variations  of  jet  flow  from  a  sonic  nozzle 


Subsonic  jet:  P^/Pgc  =  .552 


Moderately  underexpanded  jet:  p  /p  =  1.42 


Highly  underexpanded  jet:  p  /p^  =  3.57 


Figure  2.  Schlieren  photographs  of  typical  jet  types  shown  in  Figure  1. 
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2,1,2‘  Moderately  unde rexpanded  .jet.  When  the  sonic,  or 
critical,  pressure  ratio  is  reached,  a  very  weak  normal 
shock  forms  at  the  exit.  This  shock  diminishes  in  size 
rapidly  with  increasing  pressure  ratio,  however,  and  at 
V^cc  ~  1'1  the  familiar  pattern  of  "shock  diamonds"  or 
"cells"  composed  of  intersecting  oblique  shocks  is  estab¬ 


lished  in  the  core.  Except  for  a  lengthening  and  broadening 
of  the  first  few  cells  as  jet  pressure  ratio  is  increased, 
this  structure  persists  until  p/p^  *  2/  The  term  "moder¬ 
ately  underexpanded"  is  used  herein  to  denote  jets  within 
this  pressure  ratio  interval  (l.l  <  p/p^  <  2).  Because  of 
the  additional  expansion  required  in  the  unconfined  jet  flow 
beyond  the  nozzle,  the  boundaries  of  what  was  once  the 
potential  core,  in  the  subsonic  case,  are  now  determined  by 
the  requirement  of  pressure  equilibrium  between  the  outer¬ 
most  portion  of  the  flow  within  the  shock  structure  and  the 
surrounding  ambient  air.  The  initial  underexpanded  condi¬ 
tion  and  the  accompanying  shock  pattern  result  in  a  flow 
which  soon  becomes  overexpanded  at  a  point  in  the  central 
portion  of  each  cell,  in  this  region,  the  local  Mach  number 
exceeds  that  which  would  obtain  in  a  properly  expanded  jet 
with  the  same  pressure  ratio  p/p^.  The  inward  diffusion 
of  the  mixing  region,  however,  does  continue,  although  to  a 
relatively  lesser  degree  at  first,  and  ultimately  results 
in  the  complete  dissipation  of  the  shock  dominated  core. 

(In  the  absence  of  viscous  and  shock  effects,  the 
flow  would  continue  a  sequence  of  expansion  to  overexpansior 
and  recompression  to  underexpansion.)  Because  of  the  gradi¬ 
ents  of  pressure,  density,  and  Mach  number  that  exist  in 
the  core,  the  impingement  of  this  portion  of  such  a  jet 


shock  disk°rfaDnearIafoJ  P^Po°  ~  2'  at  which  the  normal 

(steyteCgUfn]f  VaSfd6d  ^^nt^Sy^eS^SthSre 


7 


might  be  expected  to  result  in  surface  pressure  distributions 
that  are  quite  sensitive  to  impingement  distance.  Downstream 
of  the  core,  of  course,  after  the  jet  has  become  subsonic, 
the  impingement  behavior  should  be  similar  to  that  of  a 
totally  subsonic  jet. 

Within  the  moderately  underexpanded  range,  effects  due 
to  instabilities  in  the  over-all  jet  flow  field  are  usually 
found  to  have  an  increased  influence  in  determining  its  decay 
and  spreading  characteristics,  both  actual  and  measured. 

This  problem  is  discussed  in  detail  in  2.2.4. 

2.1.3.  Highly  underexpanded  jet.  At  a  pressure  ratio  p-j/p^ 
of  aporoximately  2,  the  form  of  the  shock  structure  in  the 
initial  cell  begins  to  change.  Along  the  centerline,  where 
the  expansion  is  a  maximum,  the  pressure  becomes  so  low  (or 
the  Mach  number  so  high)  relative  to  ambient  pressure  that 
the  recompression  possible  in  the  remainder  of  the  cell 
through  the  existing  oblique  shocks  is  insufficient  to  raise 
the  pressure  (or  lower  the  Mach  number)  to  the  required 
initial  level  at  the  end  of  the  cell.  In  order  to  provide 
the  required  compression,  a  normal  shock  disk  forms  on  the 
centerline.  As  the  pressure  ratio  P-j/p^  is  further  increased, 
this  normal  shock  increases  both  in  strength  and  diameter. 

At  the  same  time,  the  original  oblique  shock  structure  is 
maintained  in  the  periferal  region,  although  altered  some¬ 
what  in  strength  and  shape  due  to  the  additional  expansion 
required  and  the  presence  of  the  normal  shock.  For  very 
high  pressure  ratios,  the  normal  shock  dominates  the  struc¬ 
ture  of  the  first  cell.  For  example,  with  p-i/p^  ~  20,  it 
comprises  about  40  per  cent  of  the  total  cross-sectional 
area  within  the  jet  boundaries  [ill.  It  has  also  been 
found  [ill  that  the  pressure  ratio  at  which  the  normal  shock 
disk  reappears  is  not  invariant  with  exit  Mach  number  and 
nozzle  angle.  The  value  p^/p^  ~  2  applies  only  to  a  sonic 
nozzle . 
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Immediately  downstream  of  the  normal  shock,  the  flow 
is  subsonic.  Since  the  surrounding  flow  in  the  oblique 
shock  region  remains  supersonic,  a  slip  line  exists  at  the 
boundary  between  the  two  concentric  regions.  For  a  fairly 
high  degree  of  under  expans  ion,  say  p^/p^,  ~  the  central 
subsonic  region  is  quickly  accelerated  so  that  approximately 
sonic  conditions  prevail  near  the  beginning  of  the  second 
cell.  In  this  case,  the  second  cell  may  resemble  the  first 
and  even  require  its  own  normal  shock.  For  very  high  pres¬ 
sure  ratios,  the  structure  just  downstream  of  the  first 
cell  is  not  well  defined  (for  a  recent  investigation,  see 
[17]).  However,  it  is  probable  that  tne  jet  will  be 
dominated  for  some  distance  downstream  by  the  very  strong 
normal  shock  in  the  first  cell.  Ultimately,  it  decays 
through  a  structure  with  only  oblique  shocks.  The  mixing 
■region  surrounds  the  core  as  usual,  but  its  radial  diffusion 
rate  is  small  at  first  with  the  result  that  the  effective 
core  of  the  highly  underexpanded  jet  can  be  extremely  long. 

It  should  be  noted  that  while  a  strict  definition  of  core 
length  for  any  underexpanded  jet  may  be  given  as  the  point 
at  which  the  shock  structure  disappears,  effects  due  to  the 
instabilities  present  make  this  point  difficult  to  define 
for  a  real  jet.  The  downstream  behavior  in  such  cases, 
therefore,  is  best  given  in  terms  of  the  point  at  which  the 
core's  influence  ceases.  As  in  the  case  of  other  jet 
strengths,  this  may  be  taken  as  the  point  beyond  which 
velocity  profiles  are  self  similar. 

2.2.  Experimental  program. 

Since  each  of  the  three  major  jet  variations  described 
above  was  expected  to  exhibit  an  impingement  behavior  some¬ 
what  different  from  the  others,  a  typical  case  representative 
of  each  regime  was  chosen  for  detailed  study.  Values  of 
radial  spread  and  axial  decay  for  each  of  these  jets  were 
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used  later  to  correlate  the  results  of  impingement  measure¬ 
ments  made  using  jets  with  these  same  strengths.  The  jets 
used  to  provide  these  correlation  data  are  listed  below. 

The  pressure  ratio  for  each  is  given  in  two  ways,  viz. 
p«/psc  and  pi/fV  where  Poo  ls  the  ambienfc  pressure  into 
which  the  jet  exhausts,  PgC  is  the  stagnation  pressure  in 
the  settling  chamber,  and  p1  is  the  static  pressure  (assumed 
to  be  constant)  in  the  jet  exit  plane. 


Subsonic  jet: 

-  •<&* 

pi/p«  =  1 

M]_  =  .52 


Moderately  underexpanded  jet: 
po/psc  =  *372 
Pl/P*,  =1.42 


Highly  underexpanded  Jet: 

P«/P°c  -  -148 
pl/p«>  =3.57 


Velocity  profiles  of  each  of  the  above  listed  jets  were 
measured  at  several  axial  locations.  These  locations  were 
chosen  to  represent  each  of  the  regions  of  basically  differ¬ 
ent  structure  within  a  typical  free  jet  of  high  ubsonic 
Mach  number.  They  are  listed  in  terms  of  nozzle  diameters 


downstream  from  the  exit  as  follows: 


Region  of  typical  subsonic  jet 

1.96 

Core 

7.32 

End  of  core  (transition  region) 

23.50 

Fully  developed 

39.10 

Fully  developed 

In  addition  to  the  cases  listed  above,  a  number  of 
others  were  studied  in  less  detail.  The  entire  program  is 
tabulated  in  Appendix  1. 

2.2.1.  Apparatus  and  instrumentation.  A  convergent  nozzle 
with  an  exit  diameter  d^  =  .511  inches  and  exhausting  tc 
atmospheric  pressure  was  used.  This  nozzle  was  mounted  on 
a  4.75-inch  i.d.  settling  chamber  which  was  supplied  with 
air  from  a  storage  tank  through  an  automatic  regulator 
valve.  The  maximum  storage  pressure  was  220  psig  and  the 
maximum  settling  chamber  stagnation  pressure  was  125  psig. 

The  settling  chamber  and  nozzle  are  shown  in  Figure  3,  which 
also  shows  the  flat  plate  model  used  in  a  portion  of  the 
impingement  studies  described  in  Section  3. 

Local  velocities  were  computed  from  measured  Pitot  and 
static  pressure  profiles.  The  Pitot  and  static  pressure 
probes  used  were  mounted  on  a  common  base  which  could  traverse 
the  jet  in  a  vertical  plane  at  any  axial  location  up  to  about 
60  nozzle  diameters  downstream  from  the  exit.  Both  probe 
tips  were  made  of  .032-inch  o.d.  stainless  steel  tubing;  the 
Pitot  tip  was  cut  off  square,  and  the  static  tip  was  a 
slender  ogive  with  two  .0135-inch  holes  on  opposite  sides 
5/16  inch  from  the  tip.  A  sketch  of  this  probe  and  its 
mounting  is  shown  in  Figure  4.  Pressures  were  measured  on 
liquid  manometers  or  with  Bourdon-type  test  gauges  according 
to  the  pressure  level  encountered.  Readings  for  each  run  were 
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recorded  photographically.  The  stagnation  temperature  was 
measured  with  a  bare  copper-constantan  thermocouple  in  the 
settling  chamber. 

Photographic  studies  of  the  jet  were  made  using  several 
techniques.  Schlieren  pictures  and  shadowgraphs  were  taken 
with  both  continuous  and  instantaneous  light  sources.  A 
coaxial  type  spark  source  with  a  duration  of  less  than 
1  |isec  was  used.  The  basic  optical  system  was  of  the  usual 
single  pass,  off-axis,  parallel-light  type  employing  tv, 
spherical  mirrors  of  6-inch  diameter  and  60-inch  focal 
length.  In  addition,  a  limited  number  of  pictures  were 
taken  using  a  more  sensitive  double  pass  coincident  system 
with  a  single  mirror.  While  the  resolution  of  these  latter 
pictures  is  inherently  less  than  that  achieved  with  the 
single  pass  system,  the  extra  sensitivity  provides  a  useful 
qualitative  picture  of  certain  structural  features  (see 
Subsection  2.2.3.). 

2.2.2.  Results  of  velocity  profile  and  photographic  studies. 
Results  of  the  free  jet  measurements  for  each  of  the  speci¬ 
fied  cases  chosen  for  detailed  study  are  presented  in  the 
following  paragraphs.  Basic  data  for  these  and  the  remain¬ 
ing  cases  tabulated  in  Appendix  1  are  to  be  found  in 
Appendix  2. 

For  each  typical  jet,  the  measured  total  and  static 
pressure  profiles  are  presented  with  a  spark  schlieren 
photograph  to  the  same  scale  showing  that  jet  for  the  first 
10  nozzle  diameters  downstream  (Figures  5,  7,  and  10). 

Both  total  and  static  pressures  are  plotted  in  the  form  of  a 
pressure  coefficient  expressing  the  local  value  as  a  percent¬ 
age  of  settling  chamber  gauge  pressure.  The  local  total 
pressure  is  p^  and  the  local  static  pressure  is  Pj . 
Velocity  profiles,  and  spreading  and  decay  chare.cteristics 
calculated  from  the  pressure  measurements  are  given  in 
Figures  6,  8,  and  11.  In  Figure  12,  the  spreading  character- 
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i-tics  of  all  three  Jets  are  replotted  together  in  order  to 
emphasize  certain  basic  differences. 

Local  velocities  were  computed  on  the  basis  of  the 
measured  local  pressure  ratio  and  the  measured  stagnation 
temperature  with  the  aid  of  compressible  flow  tables.  The 
pressure  ratio  at  each  point  was  evaluated  from  curves 
faired  through  the  data  for  each  pressure.*  Because  of  the 
uncertainty  in  locating  the  true  mean  axis  of  symmetry 
before  running  the  jet,  the  probes  were  traversed  through  a 
ran™  well  to  either  side  of  the  assumed  axis.  The  true 
Jet  axis  was  then  taken  to  be  the  axis  of  symmetry  of  the 
measured  profile.  The  data  were  then  replotted  with  refer¬ 
ence  to  this  true  axis  and  the  curves  drawn.  In  some  cases, 
the  radial  traverse  extended  outward  only  far  enough  to  give 
a  proper  determination  of  the  spreading  parameter  r  the 
radius  at  which  the  velocity  is  one-half  the  maximum  value. 

In  Figures  5,  7,  and  10,  a  number  of  data  points  have  been 
omitted  to  avoid  crowding.  The  velocities  shown  in  the 
profiles  of  Figures  6,  8,  and  11  are  nondimensionalized  on 
the  maximum  velocity,  even  if  it  does  not  occur  on  the 
centerline.  The  radial  coordinate  is  nondimensionalized 
on  r.5*  In  Plot  °-  decay  and  spreading  behavior,  the 
velocity  on  the  centerline  V  is  given  as  a  percentage  of 

V 

its  value  VCl  at  the  nozzle  exit.  The  centerline  value  of 
the  total  pressure  coefficient  is  also  plotted.  Specific 
structural  features  revealed  by  the  pictures  and  data  for 
each  of  the  three  jets  will  now  be  discussed. 

Subsonic  jet  (Mj  =  .52).  The  pressure  and  velocity  profiles 
of  Figures  5  and  6  clearly  reveal  the  expected  structural 
features.  The  core  with  its  profile  of  uniform  velocity  near 
the  centerline  (x/djj  =  1.96  and  3.92)  and  the  ful^y  developed 
region  with  self -similar  profiles  (x/dN  =  11.7,  23.5,39.1*  end  58.7)  have 

*For  the  subsonic  case,  the  static  pressure  was  measured 
only  for  the  three  stations  farthest  downstream.  At  other 
points,  it  was  assumed  that  p,  p  . 

,1  ~  CO 
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the  usual  appearance.*  Although  the  decay  curve  (Vc/VC]L) 
shows  the  x/d^  =  7.32  station  to  be  in  what  is  probably 
the  transition  region,  the  velocity  profile  does  not  exhibit 
any  noticeable  core  effect.  Downstream  of  x/dN  2:  11,  the 
decay  is  seen  to  follow  a  characteristic  incompressible 
l/x -dependence  quite  closely.  If  this  curve  is  extended  back 
upstream  to  a  value  of  unity,  thus  neglecting  the  transition 
region,  an  apparent  core  length  xc/d^  about  7^5  is 
found.  While  this  is  in  approximate  numerical  agreement 
with  the  results  of  other  studies  for  similar  subsonic  Mach 
numbers,  the  meaning  of  such  absolute  comparisons  is  limited, 
even  for  subsonic  .lets,  by  the  secondary  effects  already 
mentioned.  The  experiments  of  Warren  [20],  for  example, 
with  =  .69,  give  a  core  length  xc/djj  =7.2  with  a 
value  of  x/dN  ~  10  for  the  sta^t  of  the  fully  developed 
region.  Since  this  shorter  core  length  with  a  higher  Mach 
number  is  in  contradiction  with  the  usually  observed  core 
length-Mach  number  dependence,  it  is  possible  that  differ¬ 
ences  in  the  secondary  factors  affecting  the  two  experiments 
may  be  important  enough  to  account  for  this  apparent  anomaly. 
This  is  not  to  imply,  however,  that  experimental  errors 
could  not  account  for  a  difference  of  this  magnitude. 

The  spreading  behavior  is  best  observed  on  a  plot  of 
r.5/rN  aS  a  function  of  x/dN  (Figure  12).  It  is  seen 
that  the  initial  spreading  rate  decreases  slightly  for  the 
first  four  or  five  nozzle  diameters  downstream.  The  rate 
then  increases  until,  at  x/d^  ~  11*  it  becomes  fairly 
constant.  A  "transition"  region  defined  in  the  interval  in 
which  the  spreading  rate  is  changing  most  rapidly  is  seen  to 


♦Strictly  speaking,  "fully  developed"  self -similar 
velocity  profiles  imply  a  fixed  relationship  between  spread¬ 
ing  and  decay  rates.  However,  because  it  is  difficult  to 
detect  and  confirm  small  deviations  from  self-similarity  in 
the  present  data,  "fully  developed"  is  used  only  in  a 
relative  sense  in  these  discussions. 
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match  very  closely  a  region  similarly  defined  on  the  basis 
of  the  decay  curve,  i.e.  5  <  x/dN  <  11.  (  ote  that  a 

different  transition  region  is  shown  schematically  in 
Figure  1.  In  that  case  it  was  defined  as  beginning  at  the 
end  of  the  core  as  a  convenience  in  designating  the  ideal¬ 
ized  core  length.)  In  the  "fully  developed  region 
(x/dN  >  ll),  the  data  indicate  continued  slight  deviations 
from  a  truly  linear  spread.  While  it  is  felt  that  these 
deviations  exceed  experimental  error,  it  is  not  possible  to 
conclude  how  much  the  spreading  rate  actually  varies  because 
of  the  unknown  magnitude  of  jet  stability  and  turbulence 
effects.  It  is  shown  that  spreading  rates  based  on  measured 
static  pressure,  are  slightly  higher  than  those  based  on  a 
constant  ambient  static  pressure.  If  a  constant  spreading 
rate  is  determined  by  a  straight  line  fitted  among  the  three 
points  farthest  downstream  (x/d^  =  23. 5,  39.1,  and  5'3-7), 
a  spreading  angle  of  5*5°  is  found.  Using  the  ambient 
static  pressure  data,  the  angle  is  5*2°.  In  either  case, 
these  values  exceed  Warren's  result  of  4,1  by  an  amount 
that  is  probably  more  than  should  be  expected  on  the  basis 
of  the  Mach  number  difference  alone.  Warren's  data  however 
are  based  on  surveys  downstream  only  to  x/d^  =  25.  Using 
the  present  data  for  a  similar  axial  interval,  with  the 
static  pressure  assumed  equal  to  ambient  as  in  Warren's 
case,  an  angle  of  4.2°  results.  Only  if  one  can  assume 
that  differences  due  to  secondary  effects  as  well  as  Mach 
number  are  small  for  the  interval  of  axial  distance  and 
Mach  number  being  considered,  can  it  be  concluded  that  the 
agreement  is  quite  good. 

The  schli;eren  picture  shows  the  characteristic  subsonic 
turbulent  jet  mixing  region,  including  the  initial  stages  of 
the  mixing  process  just  outside  the  nozzle  exit.  The  core, 
however,  is  not  readily  discernable  because  of  the  three- 
dimensional  visual  blocking  effect  of  the  mixing  disturb¬ 
ances  (cf.  the  continuous  light  schlieren  picture  shown  in 
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Figure  2  for  the  same  case,  but  made  with  the  double  pass 
system,  in  which  the  core  is  more  easily  recognized). 

Moderately  underexpanded  Jet  ( Pi/poq=  1.42).  Effects  due  to 
under expans ion  are  at  once  apparent,  especially  in  the  axial 
decay  curve  (^ACl)  of  Figure  8.  The  centerline  velocity 
is  observed  to  be  supersonic  in  the  core  region  at  the  three 
points  chosen  for  the  measurements.  However,  because  of 
the  local  velocity  variations  to  be  expected  within  the 
length  of  each  shock  cell,  these  three  points  alone  are 
insufficient  to  show  the  detailed  core  structure,  and  the 
curve  through  them  is  thus  drawn  dashed.  Although  additional 
measurements  of  velocity  were  not  made  in  this  region,  the 
highly  detailed  survey  of  Pitot  pressure  shown  in  Figure  39 
is  indicative  of  the  kind  of  axial  variations  to  be  expected. 

The  velocity  profiles  show  clearly  the  local  effects 
of  expansion  in  the  core.  At  x/dN  =  1.96,  for  example,  the 
central  portion  of  the  profile  is  seen  to  be  supersonic.  In 
addition,  there  is  a  marked  radial  gradient  of  velocity  with 
the  peak  occurring  some  distance  from  the  centerline.  Super¬ 
sonic  central  portions  are  also  observed  for  x/dN  =  3.92 
and  ( .32,  although  the  position  of  the  peak  velocity  is 
different  in  each  case.  For  x/dN  =  11.7,  the  profile  is 
subsonic  throughout,  but  still  shows  a  slight  flattening 
near  the  centerline.  Apparently  fully  developed  subsonic 
profiles  are  observed  for  x/dN  =  23.5,  39.1,  and  58.7. 

The  behavior  of  the  measured  spreading  parameter 
r.5/r>Ni  different  from  that  observed  for  the  subsonic 
Jet  (see  Figure  12).  Beginning  at  an  axial  distance  of 
about  20  nozzle  diameters  downstream,  and  continuing  to  at 
least  40  diameters,  the  spreading  rate  in  each  axial  portion 
is  substantially  higher.  Farther  downstream,  the  rate 
decreases.  As  a  means  of  comparing  apparent  changes  in 
spreading  rate  in  different  regions  of  the  Jet,  several 
spreading  angles  have  been  computed.  In  the  interval 
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^  ^  anSle  is  only  3*0°  but  for  20  <  x/d^  <  40 

it  is  7.4  .  Using  the  two  points  farthest  downstream 
(x/dj^  =  39.1  and  58. 1),  an  angle  of  4.8°  results  for  the 
data  based  on  measured  static  pressure,  with  an  angle  of 
4.5  for  the  data  based  on  ambient  static  pressure.  It 
is  believed  that  the  increased  spreading  observed  nearest 
the  nozzle  exit  (x/d^  <  4)  is  due  to  the  widening  of  the  jet 
as  it  expands  on  leaving  the  nozzle.  Throughout  the  remain¬ 
ing  region  of  high  spreading  rate  (out'to,  say,  x/dN  =  40) 
there  is  reason  to  believe  that  the  observed  rates  are  at 
least  partly  the  result  of  jet  instability.  Photographic 
evidence  in  support  of  this  belief  is  discussed  in  2.2.3. 

The  possible  consequences  of  instability  effects  insofar 
as  the  measurements  are  concerned  can  only  be  suggested 
qualitatively  (see  2.2.4)  on  the  basis  of  the  present  data. 

Far  downstream,  the  return  to  a  lower  spreading  rate  more 
typical  of  an  incompressible  flow  appeal’s  to  be  consistent 
with  expected  trends. 

Of  particular  interest  in  the  pressure  distributions, 
shown  in  Figure  J,  is  the  behavior  of  the  static  pressure  in 
the  core  region  (also  see  Figures  AII-9  and  10).  At 
x/dN  =  a  strong  radial  gradient  is  observed,  with  a 

maximum  pressure  on  the  centerline  which  is  considerably 
higner  than  ambient,  and  a  minimum  pressure  near  the  edges 
of  the  jet  which  is  lower  than  ambient.  At  points  farther 
downstream,  the  centra:,  peak  remains,  but  the  over-all 
pressure  level  in  the  core  drops  below  ambient  pressure. 

Finally,  at  a  point  beyond  the  end  of  the  core,  the  central 
peak  disappears  and  the  over-all  level  gradually  increases 
toward  the  ambient  value.  While  this  behavior  is  qualita¬ 
tively,  both  axially  and  radially,  the  same  as  that  found  to 
exist  in  subsonic  and  properly  expam  2d  supersonic  jets 
(see,  for  example,  [.20]),  a  comparison  with  the  present  case 
can  be  misleading  without  further  clarification.  An  axial 
survey  of  centerline  static  pressure  was  made,  therefore,  to 
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help  in  understanding  this  situation.  Figure  9  presents 
the  results  of  this  survey  for  a  subsonic  jet  (p^p^,  = 

.552)  as  well  as  the  underexpanded  jet  (p,/p  =  1.42)  in 

question.  The  core  shock  structure  for  the  latter  case  is 
sketched  to  scale  so  that  the  pressures  may  be  referred  to 
their  approximate  locations  in  the  jet.  The  dashed  curve 
interpolated  among  the  data  is,  of  course,  only  a  qualita¬ 
tive  suggestion  of  the  actual  behavior.  As  such  it  is  based 
not  only  on  the  measured  points,  but  also  on  the  assumption 
that  minimum  and  maximum  pressures  occur  near  the  center  and 
end  points,  respectively,  of  each  cell.  In  any  case,  the 
extreme  gradients  within  the  core  are  clearly  evident,  and 
it  is  seen  that  the  values  found  during  profile  measurements 
(solid  symbols)  cannot  be  interpreted  as  indicating  a 
smooth  variation  in  the  axial  direction.  Also  of  interest 
is  the  fact  that  the  measured  pressure  at  the  center  of 
the  jet  exit  plane  ((p1  -  pJA FgC  “  P*,)  *  *291  at  x/dN  =  0) 
is  higher  than  that  indicated  by  the  nominal  pressure  ratio 
Pl/Poo  =  1'h2  (°r  (Pi  ‘  P«)/(Pgc  "  Pj  *  .247).  In  addi¬ 

tion,  the  axial  variation  near  the  end  of  the  core  and 
farther  downstream  closely  resembles  the  typical  subsonic 
behavior.  It  is  therefore  suggested  that  the  extreme  pres¬ 
sure  gradients  due  to  shocks  are  modified  by  a  superimposed 
radial  and  axial  distribution  which  is  similar  to  that 
existing  in  a  subsonic  or  properly  expanded  supersonic 
turbulent  jet.  Velocit"  profiles  determined  from  the  pres¬ 
sure  ratios  at  an  axial  station  in  the  shock  structure 
would  thus  reflect  a  combination  of  two  effects  (resulting, 
e.g.  in  a  peak  velocity  off  the  centerline),  and  would,  of 
course,  be  expected  to  vary  in  shape  from  point  to  point 
along  the  axis.  The  profiles  shown  in  Figure  8  for  the 
core  stations,  therefore,  are  not  necessarily  representative 
of  a  smooth  transition  of  profile  shape  from  one  axial  loca¬ 
tion  to  another. 
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The  spark  photograph  of  this  jet,  shown  In  Figure  7, 
reveals  some  distortion  of  the  stable  core  structure  as 
early  as  the  second  cell.  Farther  downstream,  the  core 
becomes  highly  unstable,  and  the  shock  cells  appear  to 
break  up  and  diffuse  Into  the  surrounding  mixing  region. 

In  the  continuous  light  picture  of  Figure  2,  which  portrays 
the  time-average  appearance,  the  downstream  cells  are  more 
easily  recognized.  Weak  sound  waves  emanating  from  the 
mixing  region  can  be  detected  in  the  spark  picture. 

Highly  underexpanded  jet  (p^p*,  =  3.57).  It  is  clear  from 
the  velocity  profile  and  decay  data  (Figure  10  and  11), 
that  the  distinguishing  structural  feature  of  this  jet  in 
the  upstream  region  is  the  normal  shock  disk  in  the  first 
cell.  The  picture  shows  that  this  shock  occurs  at  x/dN  * 
1.58.  Just  downstream  of  this  point,  at  x/d^  =  1.96,  the 
velocity  profile  exhibits  the  expected  subsonic  central 
region.  Within  this  region,  the  minimum  velocity  appears 
to  occur  just  inside  the  slip  boundary,  while  the  peak 
subsonic  velocity  lies  on  the  axis .  In  the  surrounding 
region  of  supersonic  flow,  a  peak  Mach  number  of  1.9  is 
reached,  which,  coincidentally,  happens  to  correspond  to  the 
Mach  number  for  proper  isentropic  expansion  to  p-i/p!?  =  .148. 

It  is  felt,  however,  that  at  a  point  somewhat  upstream  of 
this,  an  even  higher  Mach  number  associated  with  an  over- 
expanded  condition  should  exist.  The  photographs  "of  Figure 
10  and  Figure  2  both  reveal  an  apparent  normal  shock  in  the 
second  cell  at  x/dN  =3.3.  Slightly  downstream  of  this 
point,  at  x/dN  =  3.92,  the  velocity  profile  again  shows  a 
subsonic  central  region,  although  the  radial  extent  is  much 
less  than  it  is  for  the  x/dN  *1.96  case.  At  x/dN  =  7.32, 
the  entire  central  region  is  supersonic*  but  the  maximum 
velocity  still  does  not  occur  on  the  centerline.  In  this 
respect,  the  profile  is  similar  to  some  of  those  found  in 
the  core  of  the  moderately  underexpanded  jet.  A  substantial 
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supersonic  core  remains  at  x/djj  ~  although,  at  least 

at  this  specific  point,  the  peak  velocity  lies  on  the  axis. 
While  the  oblique  shock  structure  in  the  region  between  the 
last  normal  shock  and  the  end  of  the  core  may  not  consist 
of  well-defined  cells  such  as  those  found  in  the  moderately 
underexpanded  case,  it  is  likely  that  whatever  shocks  are 
present  will  produce  local  periodic  changes  in  the  velocity 
profile  as  long  as  they  are  of  sufficient  strength.  Thus, 
no  smooth  variation  from  profile  to  profile  should  be  inferred 
from  the  data  presented  for  this  region.  Farther  downstream, 
it  is  observed  that  the  centerline  velocity  is  just  subsonic 
at  x/dN  =  23.5.  Reference  to  the  spreading  parameter 
(r  5AN)  behavior  and  the  velocity  profiles  for  the  x/dN  = 

23 's,  39.1,  and  58.7  stations  reveals  that  a  fully  developed 
jet  flow  may  not  occur  short  of  at  least  3°  or  40  nozzle 
diameters  downstream. 

The  results  of  a  highly  detailed  Pitot  pressure  survey 
on  the  centerline  of  this  jet  are  given  in  Figure  39.  This 
survey  is  indicative  of  the  local  effects  due  to  the  normal 
shocks  present  and  the  subsequent  oblique  shock  structure  in 
the  core.  Of  particular  interest  is  the  substantial  recovery 
of  Pitot  pressure  relatively  far  downstream. 

In  order  to  verify  the  presence  of  a  normal  shock  disk 
in  the  second  cell,  some  additional  Pitot -static  pressure 
measurements  were  made  on  the  centerline  at  selected  points 
in  the  region  of  interest.  The  Mach  number  distribution 
resulting  from  these  measurements  is  shown  in  Figure  10. 

The  subsonic  region  just  downstream  of  each  normal  shock  is 
apparent.  It  is  interesting  to  note  the  sharp  increase  in 
Mach  number  from  .45  to  at  least  1.2  just  upstream  of 
the  second  shock. 

The  velocity  spread  data  for  this  jet  (Figure  12) 
reveal  a  somewhat  erratic  behavior.  In  the  region  immediate¬ 
ly  downstream  of  the  nozzle  exit,  the  bulge  observed  is 
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consistent  with  the  boundary  shape  assumed  by  the  expanding 
flow  in  the  first  few  shock  cells.  Except  for  slight  devia¬ 
tions,  the  spreading  rate  is  then  essentially  constant  for 
a  considerable  distance  downstream  (x/d^  ~  40).  In  the 
interval  4  ^  x/d^  £  12,  sPreadinS  angle  is  2.5°, 

and  for  12  <  x/dN  £  40,  an  angle  of  3*9°  is  found  using 
the  data  for  measured  static  pressure.  The  very  low  spread¬ 
ing  angle  for  4  £  x/djj  £  12  is  in  agreement  with  the  slight 
apparent  spread  observed  in  the  schlieren  picture  for  this 
case  in  Figure  2.  Downstream  of  x/d^  s:  40,  the  measured 
spread  increases.  Although  an  increase  in  this  region  seems 
to  be  consistent  with  the  appearance  of  essentially  fully 
developed  velocity  profiles  at  x/<dN  =  39.1  and  58.7,  the 
angle  of  6.0°,  based  on  the  two  data  points,  is  somewhat 
higher  than  might  be  expected  for  such  a  region  of  subsonic 
decay. 

The  continuous  light  schlieren  picture  of  this  Jet  in 
Figure  2  reveals  a  structure  downstream  of  the  second  cell 
that  seems  to  differ  somewhat  from  the  relatively  well-defined 
oblique  shock  cells  observed  at  lower  pressure  ratios. 

Although  oblique  shocks  appear  to  be  present,  the  structure 
is  more  like  that  of  a  properly  expanded  supersonic  Jet  with 
Mach  waves  in  its  core. 

2.2.3.  Special  schlieren  study  of  underexpanded  Jet.  A 
series  of  continuous  light  schlieren  pictures  was  taken 
using  the  high  sensitivity  double-pass  coincident  optical 
system.  In  this  series,  the  Jet  pressure  ratio  p^/pro  was 
varied  in  Bmall  increments  through  a  range  from  1.00  to 
more  than  4.  In  Figure  13,  a  selection  of  these  pictures 
is  shown  in  order  of  increasing  pressure  ratio.  (A  subsonic 
caBe  is  shown  for  reference.)  It  is  observed  that  an  intensi¬ 
fied  image  of  the  mixing  region  is  obtained.  Because  of  the 
relatively  long  exposure  time  (l/50  sec.),  this  image  is 
representative  of  the  time  average  appearance.  It  is  at  once 
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Figure  13,  Free  jet  from  convergent  nozzle  photographed  with  double-pass  schlieren  system. 


apparent  that  there  is  considerable  variation  of  the  observed 
spread  of  the  jet  as  the  pressure  ratio  is  changed.  This 
variation  is  interpreted  as  being  due  to  changes  in  the 
stability  characteristics  of  the  entire  jet  flow  field. 
Although  the  Instantaneous  details  of  the  structural  degrada¬ 
tion  of  the  jet  due  to  instabilities  are  readily  observed  in 
spark  pictures,  the  intensity  or  amplitude  of  the  motion  is 
difficult  to  interpret  from  single  pictures  because  of  its 
three-dimensional  nature.  The  pictures  of  Figure  13,  there¬ 
fore,  are  useful  in  making  qualitative  comparisons  of  over¬ 
all  stability  effects  for  jets  of  different  strengths.  It 
has  been  demonstrated  C 19 3  that  the  stability  of  a  given  jet 
can  depend  not  only  upon  the  pressure  ratio,  but  also  upon 
geometric  or  interference  effects  as  well  as  the  cross - 
coupling  of  acoustic  disturbances  generated  within  the  jet. 
Because  of  this  dependence,  it  is  probable  that  the  changes 
in  stability  observed  in  Figure  13  are  unique  for  this 
particular  test  apparatus.  As  an  example  of  this  uniqueness, 
it  has  been  found  that  the  proximity  of  the  mirror  (about 
4")  to  the  Jet  in  this  optical  setup  produces  a  shift  in 
what  is  believed  to  be  a  region  of  high  instability  in  the 
P/Peo  =  1.-42  jet.  A  curve  of  spreading  parameter  measured 
with  the  mirror  in  place  is  given  in  Figure  AII-14  in  order 
to  illustrate  this  effect.  The  corresponding  change  in  the 
decay  curve  (Figure  AII-16),  however,  is  relatively  small. 

This  is  consistent  with  the  assumptions  about  stability 
effects  on  profile  measurements  discussed  in  2.2.4. 

As  the  pressure  ratio  is  increased,  two  distinct  ranges 
are  noted  in  which  the  instability  appears  to  be  very  intense. 
An  increase  of  pressure  ratio  from  1.13  to  1.42  and  then 
further  to  1.5S,  spans  the  first  such  range,  with  the 
P-j/Pc  =  1.59  case  appearing  to  be  relatively  stable.  A 
second  range  of  even  greater  instability  seems  to  center 
about  the  case  for  p1/poo  =  1.84.  At  a  pressure  ratio  of 
2.00,  the  normal  shock  disk  is  first  observed.  Within  the 
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region  covered  by  the  pictures,  the  degree  of  instability 
appears  to  lessen  with  further  increases  in  pressure  ratio 
above  2.  It  should  be  noted  that  although  the  Jet  chosen 
(P-/P*,  =  1.42)  for  detailed  study  in  the  moderately  under¬ 
expanded  case  appears  to  fall  within  one  of  the  ranges  of 
high  instability,  the  behavior  shown  in  the  picture  is  not 
in  itself  conclusive  because  of  the  aforementioned  mirror 
proximity  effect. 

The  time -average  appearance  of  the  core  shock  structure 
is  also  of  interest.  As  might  be  expected,  the  more  unstable 
Jets  show  fewer  well-defined  cells.  The  contrast  is  particu¬ 
larly  great  between  the  cases  for  p^/p^  -  1.59  and  1.84. 
Once  the  normal  shock  occurs  in  the  first  cell,  there  is  a 
gradual  change  in  the  appearance  of  cells  farther  downstream. 
At  first,  these  cells  seem  to  follow  the  characteristics -type 
of  pattern  used  as  a  model  for  the  moderately  underexpanded 
Jet.  However,  between  the  pressure  ratios  2.59  and  3*57, 
the  regular  cellular  division  seems  to  give  way  to  a  more 
continuous  pattern  of  Intersecting  oblique  shocks. 

2.2.4.  Factors  affecting  profile  measurements.  The  velocity 
profiles  upon  which  the  Jet  spreading  and  decay  results  are 
based  were,  of  course,  determined  from  measurements  made  with 
Pitot  and  static  pressure  probes.  Inherent  in  such  measure¬ 
ments  are  certain  limitations  introduced  by  the  properties 
of  the  flow  itself.  In  the  case  of  turbulent  Jets,  the  most 
important  limiting  factors  are  thought  to  be  the  turbulence 
in  the  mixing  region  and  the  over-all  Jet  instability,  (it 
is  felt  that  alignment  errors  due  to  neglect  of  the  radial 
component  of  the  mean  velocity  and  probe  "angle -of -at tack 
errors  due  to  the  shear  flow  mean  profile  are  of  minor 
importance.)  Although  no  quantitative  evaluation  of  these 
factors  is  possible  on  the  basis  of  the  existing  data,  some 
general  conclusions  about  the  relative  validity  of  the 
measurements  should  be  possible. 


The  pressure  sensed  by  a  static  pressure  orifice  is 
affected  by  transverse  velocity  components  arising  from 
turbulence  as  well  as  any  other  phenomenon  having  a  cross¬ 
wise  component.  The  magnitude  and  sign  of  the  resulting 
error,  however,  depend  on  a  complex  relationship  among 
probe  size,  turbulence  scale,  and  the  magnitude  and  space 
correlation  of  local  velocity  fluctuations.  In  general, 
therefore,  the  validity  of  the  static  pressure  measurements 
can  be  assumed  to  be  the  greatest  in  regions  of  the  Jet 
where  turbulence  and  instability  are  the  least  relative  to 
the  magnitude  of  the  mean  motion,  namely,  in  the  upstream 
regions . 

In  many  cases,  jet  velocity  profiles  are  determined 
from  the  measurement  of  Pitot  pressures  alone,  with  the 
static  pressure  considered  to  be  constant  and  equal  to 
ambient  pressure.  In  the  present  experiments.  In  which 
static  pressures  were  measured  in  most  cases, it  is  possible 
to  compare  profile  parameters  determined  in  both  ways.  In 
Figure  12,  values  of  the  spreading  parameter  based  on  ambient 
static  pressure  are  shown  for  several  cases.  It  is  seen 
that  a  somewhat  larger  spreading  rate  results  when  measured 
static  pressures  are  used.  It  is  not  possible,  however,  to 
determine  the  degree  to  which  these  measured  pressures 
actually  contribute  to  the  determination  of  a  true  profile 
because  the  measurements  are  most  in  doubt^wirere^  -tfiSy  can 
have  the  greatest  influence,  i.e.  in  the  outer  portion  of 
the  downstream  region  where  they  approach  the  magnitude  of 
the  total  pressure  level. 

While  total  or  Pitot  pressure  measurements  are  also 
affected  by  turbulent  velocity  components,  it  is  felt  that 
over-all  Jet  instability  effects  may  be  of  greater  impor¬ 
tance  in  some  of  the  present  cases.  In  such  cases,  the 
response  of  the  Pitot  tube  at  each  point  in  the  profile  can 
be  considered  to  be  that  resulting  from  a  fluctuating  veloc¬ 
ity  at  that  point.  If  It  is  assumed  that  this  response 


represents  the  time-average  value  of  the  fluctuation,  a 
typical  Jet  mixing  profile  measured  in  this  way  will  differ 
from  its  instantaneous  shape.  Assuming  a  lateral  disturb¬ 
ance  motion  whose  mean  amplitude  is  distributed  axisymmet- 
rically,  the  measured  profile  would  appear  to  be  somewhat 
flattened  at  the  center  and  spread  out  at  the  edges. 

Spreading  and  decay  rates  based  on  such  time-average  profiles 
would,  of  course,  be  larger  than  those  based  on  instantane¬ 
ous  profiles. 

Because  of  the  foregoing  factors,  it  is  clear  that  the 
measurement,  for  example,  of  a  high  spreading  rate  for  a 
given  Jet  may  only  be  indicative  of  the  fact  that  the  Jet  is 
highly  unstable.  The  instability  would  then  have  to  be 
either  eliminated  or  evaluated  by  other  means  before  the 
true  viscous  spreading  rate  could  be  determined. 

2.3.  Discussion  of  results  and  comparison  with  theory. 

The  main  objective  of  the  foregoing  study  of  free  Jet 
properties  has  been  the  determination  of  spreading  and  decay 
characteristics  to  be  used  to  correlate  the  results  of 
impingement  experiments  using  the  same  Jet  apparatus.  By 
using  such  data  in  this  way,  the  influence  on  the  correla¬ 
tion  of  secondary  effects  such  as  Jet  stability  might  be 
expected  to  be  minimized.  Also  of  interest  in  this  study 
has  been  the  general  behavior  and  structure  of  free  Jets 
themselves,  especially  cases  in  which  the  Jet  is  under¬ 
expanded  . 

The  results  of  the  three  typical  cases  presented  in 
detail  in  2.2.2  confirm  a  number  of  expected  similarities  as 
well  as  important  differences  among  the  basic  flow  types. 

In  the  core  region  of  each  Jet,  the  differences  are  most  in 
evidence.  The  core  of  the  subsonic  Jet  is,  of  course, 
determined  by  the  inward  diffusion  of  the  turbulent  mixing 
region,  whereas  the  moderately  underexpanded  Jet  has  an 
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additional  determining  influence  in  the  system  of  oblique 
shocks  present.  For  the  highly  underexpanded  case,  the 
normal  shock  disk  is  a  dominant  factor  in  the  local  struc¬ 
ture  of  the  core.  Because  of  the  very  presence  of  shocks 
in  the  underexpanded  cases,  however,  it  is  difficult  to 
specify  a  consistent  criterion  for  core  length  that  can  be 
applied  with  equal  pertinence  to  all  the  jets.  It  is  felt, 
therefore,  that  the  most  meaningful  basis  of  comparison  is 
the  downstream  behavior  in  terms  of  the  point  at  which  a 
fully  developed  turbulent  mixing  profile  is  observed.  The 
present  data  are  sufficiently  detailed  to  be  used  in  this 
way. 

Using  the  measured  velocity  profiles  by  themselves,  it 
is  found  that  the  subsonic  jet  can  be  considered  fully 
developed  somewhere  between  x/dN  =  7.32  and  11.7.  The 
constant  relationship  between  centerline  velocity  decay  and 
jet  width  or  spread,  which  is  thus  implicit  and  which  must 
hold  if  axial  momentum  is  to  be  conserved,  is  fairly  well 
confirmed  in  separate  plots  of  these  two  parameters.  The 
moderately  underexpanded  jet  (p^/p^  =  1.42),  however, 
exhibits  a  profile  at  x/dN  =  11.7  that  still  does  not 
match  those  far  downstream.  It  has  been  pointed  out  that 
this  particular  jet  appears  to  be  quite  unstable  and  that 
measured  velocity  profiles  may  represent  a  distortion  of 
the  actual  instantaneous  profile.  Because  of  this,  the 
definition  of  a  fully  developed  region  is  difficult.  It 
is  seen,  for  example,  that  the  velocity  profiles  are  very 
close  to  being  similar  for  x/dN  =  23.5,  39.1,  and  58.7, 
while  at  the  same  time  there  is  a  marked  decrease  in  the 
spreading  rate  in  the  same  range.  This  situation  could 
result  if  jet  instabilities  were  stronger  in  the  upstream 
region  and  thus  resulted  in  broader  measured  profiles  there. 
In  fact,  if  it  is  assumed  that  the  measurements  far  down¬ 
stream  at  x/dN  =  58.7  are  relatively  unaffected  by  instabi¬ 
lity,  it  is  found  that  the  over-all  spreading  rate  required 
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to  reach  the  measured  width  at  that  point  is  very  nearly 
the  same  as  that  required  for  the  subsonic  jet  at  the  same 
point.  The  highly  underexpanded  jet  is  apparently  dominated 
by  a  very  long  supersonic  core,  as  shown  by  the  low  initial 
spreading  rate  and  the  centerline  Mach  number  survey.  It  is 
doubtful,  in  fact,  that  a  fully  developed  region  occurs  at 
all  within  the  range  of  the  measurements.  Velocity  profiles 
for  x/d^  =39.1  and  58.7  are  essentially  similar,  but  the 
spreading  rate  between  these  points  is  higher  (6.0°)  than 
that  usually  associated  with  a  fully  developed  subsonic 
mixing  region. 

A  comparison  has  been  made  between  the  results  of  this 
study  and  the  semi -empirical  integral  analysis  of  Warren  [20] 
based  on  Prandtl's  constant  exchange  coefficient  concept. 

This  theory  differs  from  the  usual  mixing  length  hypothesis 
in  that  it  defines  the  effective  eddy  viscosity  or  exchange 
coefficient  €  directly  in  terms  of  the  mean  flow  proper¬ 
ties.  For  a  typical  fully  developed  jet  mixing  region,  it 
is  assumed  that 


where  K  is  a  proportionality  constant  to  be  determined 
experimentally,  Warren  found  that  K  could  be  correlated 
with  M^j^  within  his  experimental  range.  This  correlation, 
which  was  based  on  Warren's  data  for  both  subsonic  and 
ideally  expanded  supersonic  jets,  is  given  by 

X  =  .0430  -  .0069  M2 

The  principal  objective  of  the  present  comparison  is  to  see 
how  well  the  decay  behavior  of  an  underexpanded  Jet  may  be 
correlated  with  that  for  a  Jet  that  is  properly  expanded  at 
the  same  pressure  ratio.  Although  it  is  probable  that  the 
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relationship  of  K  and  ir  unique  for  a  given  test 

condition,  the  determination  of  such  a  relationship  is  not 
within  the  scope  of  the  present  experiments.*  Therefore  the 
comparison  with  Warren's  method  is  carried  out  on  the  basis 
of  his  correlation  of  K  and  M^. 

In  order  to  compare  the  measured  decay  of  a  given  under - 

expanded  jet  and  that  computed  for  a  properly  expanded  jet  of 

the  same  pressure  ratio  p ,  it  is  first  assumed  that  the 

Jet  exit  locations  coincide.  The  nondimens ional  axial 

coordinate  x/djj  of  the  computed  jet  is  then  based  on  a 

diameter  given  by  the  area  ratio  for  proper  isentropic 

expansion  to  the  given  pressure  ratio  with  the  throat  area 

equal  to  that  of  the  actual  nozzle.  Similarly,  the  ratio  of 

the  exit  velocity  of  the  equivalent  properly  expanded  jet  to 

that  measured  at  the  sonic  exit  of  the  underexpanded  jet  is 

used  to  scale  the  entire  computed  decay  curve.  The  results 

of  this  type  of  comparison  for  several  cases  are  given  in 

Figure  14.  It  is  seen  that  the  degree  of  correlation  is  very 

good  for  the  subsonic  Jet,  but  somewhat  varied  for  the  other 

cases.  In  Figure  15,  the  ratio  of  measured  to  computed  decay 

parameter  (V  /V  )  is  plotted  as  a  function  of  jet 
Gexp  Gth 

pressure  ratio  in  the  underexpanded  regime  for  several  axial 
locations.**  Within  a  range  of  pressure  ratios  centered 
about  that  for  the  formation  of  the  normal  shock,  the  agreement 
is  no  better  than  60-70  per  cent.  For  pressure  ratios 
p^/pw  >  3,  however,  the  agreement  is  much  better.  It  is  also 

*In  fact,  it  is  possible  that  a  better  correlation  can  be 
found  if  K  is  assumed  to  be  a  function  of  some  local  Mach 
number  which  is  characteristic  of  the  flow  at  each  axial  station. 
Three  such  Mach  numbers  which  have  been  suggested  are  those  on 
the  centerline,  on  the  dividing  streamline,  and  on  the  stream¬ 
line  at  r  c. 

•  o 

**The  pressure  ratio  for  appearance  of  the  normal  shock 
shown  in  Figure  15  was  determined  from  a  plot  of  shock  diameter 
as  a  function  of  pressure  ratio  by  extrapolating  the  shock 
diameter  to  zero. 


that  given  by  application  of  Warren 
pressure  ratio  and  a  Mach  number  M 


Figure  15.  Ratio  of  measured  centerline  velocity  to  that  found  for  equivalent 

jet  pressure  ratio. 
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observed  that  while  the  per  cent  agreement  falls  off  with 
axial  distance,  the  difference  in  agreement  due  to  the  pressure 
ratio  effect  is  also  less  far  downstream.  Since  the  theory- 
does  not  account  for  secondary  effects,  these  results  lead  to 
conclusions  that  are  quite  consistent  with  those  based  on 
photographic  evidence  and  measured  spreading  rates  alone, 
namely,  that  jet  instability  effects  may  be  large  in  the 
moderately  underexpanded  range,  and  that  such  effects  are 
diminished  at  points  far  downstream.  It  is  also  evident  that, 
except  in  the  moderately  underexpanded  range,  the  semi- 
empirical  method  used  for  a  properly  expanded  jet  results  in 
a  reasonably  good  approximation  to  the  underexpanded  jet 
within  the  range  of  pressure  ratios  investigated.  While  the 
core  shock  structure  can,  tinder  certain  conditions,  have  a 
definite  influence  on  the  stability  of  the  jet,  it  is  apparent 
that  its  over-all  effect  on  decay  rates  is  minor. 

It  should  be  pointed  out  that  the  known  applications  of 
Warren's  method  to  properly  expanded  jets  have  usually  been 
restricted  to  axial  distances  of  less  than  x/dN  =  25.  For 
this  reason,  it  is  not  possible  to  verify  that  the  method  is 
any  better  for  properly  expanded  jets  far  downstream  than  it 
is  for  the  present  underexpanded  ones.  The  single  point  for 
=  .515  at  x/dN  =39.1  shown  in  Figure  14  is,  of  course, 
by  itself  inconclusive. 
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3 .  IMPINGEMENT  STUDIES 

3.1.  Basic  flow  characteristics. 

Recent  interest  in  ground  effects  machines,  V/STOL 
aircraft,  and  the  vertical  launching  and  landing  of  rockets 
has  led  to  a  number  of  studies  of  various  aspects  of  the  jet 
impingement  problem.  In  addition,  there  has  been  increased 
study  of  certain  industrial  processes  involving  heating  by 
impinging  hot  Jets  and  flames.  In  the  ground  effects  and 
V/STOL  field,  the  need  to  understand  impingement  processes 
has  arisen  not  only  with  regard  to  increasing  the  vehicle's 
lifting  effectiveness  while  in  ground  proximity  [27-34],  but 
also  in  connection  with  downwash  erosion  effects  on  the 
ground  below  [ 35] •  Study  of  the  ground  erosion  effect  has 
also  been  extended  to  the  problem  of  landing  rockets  on 
hypothetical  lunar  and  planetary  surfaces  [36-40],  Problems 
associated  with  the  impingement  and  deflection  of  rocket 
exhausts  and  the  resulting  loading  of  adjacent  surfaces  have 
been  treated  both  theoretically  and  experimentally  [41-45]. 
The  basic  problem  of  determining  heat  transfer  between 
surfaces  and  impinging  Jet  flows  has  also  been  investigated 
in  a  variety  of  ways  [46-59].  Other  investigators  have 
emphasized  the  basic  aspects  of  flow  processes  involved  in 
impingement  [60-68]  as  well  as  certain  special  problems  such 
as  noise  generation  [69] . 

The  flow  field  produced  when  an  axially  symmetric  air 
Jet  impinges  on  a  solid  surface  held  normal  to  it  consists  of 
three  general  regimes.  First,  there  is  the  jet  Itself,  up¬ 
stream  of  the  point  where  any  local  influences  due  to  the 
strong  interaction  of  the  impingement  are  felt.  Throughout 
this  regime,  of  course,  secondary  effects  (such  as  those 
described  in  Section  2)  produced  by  the  impingement  surface 
will  undoubtedly  play  a  part  in  determining  the  exact  jet 
characteristics.  The  second  regime  of  interest  is  the 
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impingement  regime,  wherein  the  flow  properties  are  primarily 
determined  by  the  direct  interaction  of  the  jet  and  the  solid 
surface.  It  is  here  that  the  large  gradients  of  pressure, 
density,  and  velocity  associated  with  the  rapidly  changing 
flow  direction  occur.  Once  the  flow  has  been  completely 
turned  in  a  direction  parallel  to  the  impingement  surface  and 
is  no  longer  influenced  locally  by  impingement  processes,  it 
enters  the  third  basic  regime,  that  of  the  wall  jet.  Here 
the  primarily  radial  flow  develops  into  a  fully  developed 
wall  jet  characterized  by  an  inner  layer  of  boundary  layer¬ 
like  flow  and  an  outer  layer  of  free  shear  turbulent  mixing. 
It  is  probable,  of  course,  that  the  character  of  at  least  the 
first  two  of  these  regimes  will  be  highly  sensitive  to  local 
changes  in  the  structure  of  the  impinging  jet,  especially  for 
cases  in  which  an  underexpanded  jet  impinges  at  close  range. 
Each  of  the  basic  regimes  and  the  symbols  used  to  designate 
certain  quantities  are  shown  in  Figure  16. 

In  order  to  estimate  the  heat  transfer  at  the  stagnation 
point  of  an  impinging  flow  such  as  that  just  described,  the 
local  radial  velocity  gradient  (due/dr)r=0  which  appears  as 
a  parameter  in  the  usual  stagnation  point  heat  transfer 
equation  must  be  evaluated.  This  can  be  done  experimentally 
by  relating  the  parameter  to  the  static  pressure  distribution 
on  the  surface  in  the  immediate  vicinity  of  the  stagnation 
point.  Assuming  the  flow  outside  the  boundary  layer  to  be 
locally  incompressible,  the  local  pressure  in  the  laminar 
boundary  layer  may  be  written 

P  -  P?  '  |p eue  (1) 


or 


dp 

dr 


dUe  1..2  f£e 

_peue  dr  2  e  dr 


(2) 


where  Pp  is  the  total  pressure  at  the  stagnation  point 
(i.e.  where  u  »  0),  ue  is  the  velocity  at  the  edge  of 
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Figure  16 .  Schematic  view  of  impingement  flow  with  flat  plate  model 

in  place . 


Cylindrical  cup 


Figure  17.  Basic  impingement  model  shapes. 
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the  boundary  layer,  and  pg  Is  the  local  density.  Neglecting 
the  second  term  on  the  right,  since  dpg/dr  -  0, 


At  r  =  0,  ue  =  0,  so  that 


Thus  the  parameter  in  question  is  proportional  to  the  square 
root  of  the  curvature  of  the  pressure  distribution  at  the 
stagnation  point.  For  the  purpose  of  evaluating  (du  /dr)  _ 
directly  from  measured  pressure  data,  however,  an  alternate 
form  of  this  relation  is  derived  directly  from  Equation  1 
making  use  of  the  equation  of  state  Pp  =  p  RT°.  Thus,  expand- 
ing  ue  in  a  power  series  about  the  stagnation  point  r  =  0, 


where  r  is  the  wetted  radius  of  the  impingement  surface, 

0  w 

T  is  the  stagnation  temperature  of  the  flow,  and  R  is  the 
specific  gas  constant.  Solving  for  (dUg/dr)r_Q 
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we  have 


The  evaluations  of  this  parameter  for  a  number  of  impinge¬ 
ment  conditions  are  given  in  3*2,3.  The  experimental  program 
is  described  in  detail  in  the  next  subsection. 

3.2.  Experimental  program. 

The  bulk  of  the  experimental  program  was  devoted  to  a 
study  of  the  normal  impingement  (a  =  90°)  of  the  three  basic 
Jet  types  described  in  Section  2.  Each  of  these  Jets  was 
impinged  upon  four  different  model  shapes  and  the  stagnation 
region  pressure  distributions  determined.  Stagnation  point 
radial  velocity  gradients  were  then  computed.  In  addition, 
two  methods  of  nondimensionalizing  the  measured  gradients  in 
terms  of  known  Jet  properties  were  evaluated.  One  such 
method  was  based  on  conditions  at  the  Jet  nozzle  exit,  and 
the  other  on  local  conditions  at  the  impingement  station  in 
the  free  Jet.  These  nondimensional  forms  are  discussed  in 
3.2.3. 

The  three  Jets  used  for  the  impingement  studies  were  as 
listed  in  the  discussion  of  the  free  Jet  program  in  2.2, 
except  for  the  subsonic  Jet  which  had  a  negligibly  different 
pressure  ratio,  i.e,  P^/p^q  =  .800  rather  than  .834.  The 
moderately  underexpanded  Jet  had  p^/p^  =  .372,  and  the 
highly  underexpanded  Jet  had  p^p^  =  .148.  The  impingement 
distances  chosen  were  also  the  same  as  those  for  the  free 
Jet  experiments,  i.e.  x/dN  =  1.96,  7.32,  23.5,  and  39.1. 

In  addition,  several  other  locations  were  used  in  order  to 
fill  in  data  in  regions  of  special  interest .  The  entire 
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program  of  normal  impingement  cases,  which  included  measure¬ 
ments  of  over-all  surface  pressure  distributions  and  certain 
additional  studies  with  the  flat  plate  model,  as  well  as  the 
detailed  stagnation  region  measurements  is  tabulated  in 
Appendix  III.  An  additional  program  devoted  primarily  to 
the  study  of  impingement  on  the  flat  plate  at  oblique  angles 
(a  <  90°)  is  to  be  reported  separately. 

3.2.1.  Apparatus  and  instrumentation.  Except  for  the  impinge¬ 
ment  models  themselves,  the  test  setup  was  exactly  the  same 
as  that  used  for  the  free  jet  studies  (2.2.1).  In  Figure  3, 
the  nozzle  is  shown  with  the  flat  plate  model  mounted  in 
position.  The  mounting  was  designed  so  that  fine  centering 
adjustment  could  be  made  either  horizontally  or  vertically. 

The  impingement  angle  was  adjusted  by  means  of  a  jack  screw 
which  rotated  each  model  about  a  horizontal  line  passing 
through  its  stagnation  point.  Axial  changes  in  impingement 
distance  were  made  by  shifting  the  entire  model  supporting 
structure  to  the  desired  location  along  two  steel  angle 
rails  at  the  bottom.  The  mounting  as  a  whole  was  made  to 
be  rigid  enough  to  minimize  deflection  under  jet  dynamic 
loading,  while  at  the  same  time  having  the  main  members  as 
far  removed  as  feasible  from  the  Impingement  region  so  as 
to  minimize  the  possibility  of  Interference  with  the  flow. 

All  of  the  impingement  models  were  made  with  the  same 
wetted  diameter,  I.e.  the  distance  along  the  impingement 
surface  from  edge  to  edge  through  the  center.  This  distance, 
based  on  a  hemisphere  model  diameter  of  6  inches  was  9.4-2 
inches.  The  individual  model  characteristics  were  as 
follows  (see  Figure  17): 

Flat  plate:  Aluminum  disk,  1/2  inch  thick  and  9.42 

inches  in  diameter.  27  pressure  taps  along  vertical 

diameter  and  15  along  horizontal  diameter. 
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Convex  and  concave  hemispheres:  Fiber  glass -epoxy- 
resin  molded  shells,  about  3/8  inch  thick  and  6  inches 
in  diameter.  Concave  model  made  first  on  male  mold; 
concave  model  then  used  as  mold  for  convex  model. 

27  pressure  taps  along  vertical  diameter  and  5  along 
horizontal  diameter. 

Cylindrical  cup:  Brass  flat  plate  1/2  inch  thick  and 
6  inches  in  diameter  with  brass  cylindrical  rim 
1/8  inch  thick  and  1.71  inches  high.  28  pressure  taps 
along  vertical  diameter  including  rim  and  11  along 
horizontal  diameter. 

All  the  pressure  tap  holes  were  drilled  with  a  number  75 
drill  except  for  four  closely  spaced  (l/l6  inch  spacing) 
holes  including  the  stagnation  point  hole  which  were  number 
80  (.0135  inch).  The  latter  holes  were  along  the  vertical 
diameter. 

Pressure  readings  for  each  run  were  made  on  multiple 
manometers  or  test  gauges  and  were  recorded  photographically. 
The  stagnation  temperature  was  measured  as  before. 

3.2.2.  Results  of  pressure  distribution  measurements  and 
photographic  studies . 

Pressure  distributions.  For  each  of  the  basic  combinations 
of  Jet  strength  and  impingement  distance,  both  detailed 
stagnation  region  and  over-all  pressure  distributions  were 
measured  for  each  model.  Of  particular  importance  in 
obtaining  these  results  was  the  initial  alignment  of  the 
models  relative  to  the  Jet  flow,  A  model  was  first  aligned 
parallel  to  the  jet  exit  plane  by  means  of  direct  measure¬ 
ment  between  its  outer  edges  and  a  straightedge  held  across 
the  nozzle  exit.  Vertical  and  horizontal  centering  were  then 
accomplished  with  the  Jet  running  by  nulling  the  pressure 
differential  between  pressure  taps  equidistant  from  the 
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center  tap.  Although  the  initial  parallel  alignment  cannot 
take  into  account  the  fact  that  the  jet  may  not  issue  from 
the  nozzle  exactly  coaxial  with  it  or  that  interference  or 
buoyant  effects  may  change  the  jet's  direction  slightly,  it 
was  found  that  stagnation  region  distributions  were  quite 
insensitive  to  small  changes  in  impingement  angle.  Thus  it 
was  felt  that  these  factors  could  be  neglected  as  long  as 
the  model  was  centered. 

In  order  to  obtain  high  resolution  distributions  in  the 
vicinity  of  the  stagnation  point,  a  technique  was  used  of 
translating  the  model  slightly  between  each  of  several  jet 
runs.  In  this  case,  the  model  was  exactly  centered  only 
horizontally.  The  vertical  centering  adjustment  was  then 
used  to  shift  the  model  slightly  up  or  down  in  its  own 
surface  plane  so  as  to  bring  the  central  pressure  taps  to  a 
different  point  in  the  flow.  Because  of  the  usually  consist¬ 
ent  data  obtained,  it  was  assumed  that  the  maximum  over-all 
translation  of  about  l/8  inch  had  little  or  no  effect  on  the 
absolute  position  of  the  jet.  A  typical  pressure  distribu¬ 
tion  resulting  from  the  use  of  this  technique  is  shown  in 
Figure  18. 

The  results  of  these  detailed  stagnation  point  measure¬ 
ments  are  shown  for  all  the  basic  combinations  in  Figures  19 
through  30.  On  these  plots  the  pressure  is  given  as  a  ratio 
of  local  to  stagnation  point  absolute  value,  while  the  radial 
distance  is  given  in  nozzle  radii  from  the  stagnation  point. 
In  order  to  show  clearly  the  slight  differences  between  the 
distributions  for  different  shapes  on  each  plot,  the  data 
points  have  been  omitted.  For  reference,  the  Pitot  pressure 
profile  of  the  free  jet  at  each  axial  location  is  also 
plotted. 

Except  for  the  case  P^/p^  =  3.57’j  x/dN  =  is 

observed  that  the  pressure  distributions  follow  the  general 
local  character  of  the  free  jet  (cf.  Figures  5,  7,  and  10), 
although  there  is  a  tendency  for  the  impingement  distribution 
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Comparison  of  measured  stagnation  region  pressure  distributions 
for  several  surface  shapes. 


Comparison  of  measured  stagnation  region  pressure  distributions 
for  several  surface  shapes . 
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Figure  24.  Comparison  of  measured  stagnation  region  pressure  distribut 

for  several  surface  shapes. 


Figure  26.  Comparison  of  measured  stagnation  region 


Figure  27.  Comparison  of  measured  stagnation  region  pressure  distributions 
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to  be  relatively  spread  out.  The  excepted  case  exhibits  a 
flattened  distribution  with  no  distinct  peaks.  Other  cases 
involving  the  highly  underexpanded  jet  show  apparent  stagna¬ 
tion  points  to  either  side  of  the  center  with  a  relatively 
flat  distribution  or  a  third  peak  at  the  center.  This  behav¬ 
ior  has  been  investigated  further  and  is  discussed  in  3.2.4. 

A  consistent  trend  is  that  due  to  surface  shape.  The  convex 
hemisphere  distributions  always  drop  off  more  rapidly  than 
those  of  the  flat  plate,  while  those  for  the  concave  shapes 
drop  off  less  rapidly.  The  relative  behavior  of  the 
cylindrical  cup  and  concave  hemisphere  shapes  appears  to 
depend  on  the  impingement  distance. 

Over-all  pressure  distributions,  measured  from  edge  to 
edge  across  the  vertical  wetted  diameter  of  each  model,  are 
shown  in  Figures  31  through  38  for  each  jet  strength  and 
model  shape.  These  data  are  plotted  with  the  local  pressure 
as  a  percentage  of  its  value  Pp  at  the  stagnation  point. 
Thus  it  is  possible  to  note  the  relative  change  in  shape  of 
the  distributions  for  different  impingement  distances.  Plots 
of  these  same  data  given  in  terms  of  the  jet  settling  chamber 
stagnation  pressure  are  included  in  Appendix  IV. 

In  general,  the  effect  of  distance  as  shown  in  Figures 
31  through  38,  is  as  might  be  expected,  with  the  pressure 
distribution  following  the  basic  spreading  trends  of  the 
free  jet.  Certain  details,  however,  are  apparent  which  are 
not  observed  in  the  free  jet  data.  In  particular,  there  is 
often  a  distinct  reversal  of  the  radial  pressure  gradient 
in  a  short  interval  between  2  and  3  nozzle  radii  from  the 
stagnation  point  when  x/d^  =  1.96  and  7.32.  This  reversal 
depends  to  some  extent  on  the  surface  shape,  being  strongest 


♦This  reversal  is  to  be  distinguished  from  that  found 
near  the  center  for  the  Pi/Poo  =  3-57  jet  which  does  occur 
in  the  free  jet  and  which  is  apparently  associated  with  the 
jet  normal  shock  disk. 
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for  the  convex  hemisphere.  Since  it  occurs  for  this  shape 
even  for  subsonic  impingement,  it  appears  not  to  be  associa¬ 
ted  specifically  with  jet  underexpansion  phenomena.  Whether 
or  not  the  reversed  gradient  is  sufficient  to  cause  local 
separation  is  impossible  to  determine  from  these  data  alone. 
Another  region  of  reversed  gradient  occurs  near  the  outer 
edges  for  all  the  surfaces  but  the  flat  plate,  although, 
again,  separation  cannot  be  confirmed. 

The  behavior  of  the  concave  shapes  (hemisphere  and  cup) 
for  impingement  far  downstream  is  of  particular  interest. 

In  these  cases,  it  is  seen  that  a  large  portion  of  the  surface 
is  subjected  to  a  pressure  nearly  as  high  as  that  at  the 
stagnation  point.  In  effect,  the  entire  flow  inside  these 
shapes  approaches  a  stagnation  condition  when  the  jet  has 
spread  to  a  size  comparable  to  that  of  the  model.  In  the 
cylindrical  cup,  the  stagnation  regions  in  the  sharp  corner 
are  quite  clear.  The  over-all  stagnating  effect  is  approached 
the  closest  for  this  model,  especially  in  the  case  of  impinge¬ 
ment  by  the  moderately  underexpanded  jet  (p^Poo  "  1.^2).  It 
is  seen  that  the  corner  stagnation  pressures  are  as  high  as 
80  per  cent  of  those  at  the  central  stagnation  points. 

The  results  of  a  highly  detailed  axial  survey  of  the 
pressure  at  the  center  of  the  flat  plate  model  are  given  in 
Figure  39.  Also  plotted  is  a  similar  survey  of  the  Pitot 
pressure  at  the  center  of  each  free  jet.  Pitot  pressure  is 
plotted  because  the  total  pressure  loss  due  to  the  stand-off 
normal  shock  cannot  be  determined  for  the  impingement  case. 
The  data  points  for  the  free  jet  have  been  omitted  for 
clarity,  a  measurement  having  been  made  at  intervals  of 
0.2  nozzle  diameters  through  most  of  the  range  0  £  x/du  <  10 • 
The  values  of  xc/dN  shown  represent  core  lengths  for  an 
equivalent  properly  expanded  jet  of  the  same  pressure  ratio 
computed  using  Warren’s  method  (see  2.3) •  For  the  subsonic 
case,  as  might  be  expected,  the  curves  are  very  close.  In 
fact,  as  plotted  on  the  present  scale,  they  cannot  be 
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Figure  39.  Comparison  of  pressures 
at  center  of  flat  plate  model 
with  Pitot  pressures  on 
free  jet  center  line. 
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separated.  Both  of  the  underexpanded  cases  reveal  the 
degree  to  which  small  local  changes  due  to  shock  structure 
are  duplicated  in  the  two  kinds  of  measurement .  The  most 
important  difference  observed  between  the  free  jet  and 
impingement  data,  is  an  axial  shift  such  that  the  plate  is 
seen  to  experience  a  pressure  which  actually  occurs  in  the 
free  Jet  at  a  point  upstream  of  the  plate  location.  This 
shifting  effect  is  thought  to  be  at  least  partly  the  result 
of  the  different  stand-off  distances  of  the  impingement 
normal  shocks .  Since  for  a  given  Mach  number  the  stand-off 
distance  of  a  normal  shock  increases  with  the  effective 
bluntness  of  the  body*  it  is  clear  that  the  probe  stand-off 
distance  must  be  considerably  less  than  that  for  the  plate. 
Other  factors  contributing  to  the  shift  may  involve  possible 
changes  in  the  stability  of  the  jet  itself  due  to  the 
secondary  interference  effect  of  the  plate  on  the  jet  flow 
field.  In  this  regard,  it  is  interesting  to  note  that  in 
the  case  for  p-^/p^  =  1.42  the  shift  continues  in  the  same 
direction  to  points  well  downstream  of  the  supersonic  portion 
of  the  jet.  Some  evidence  has  been  given  that  seems  to  show 
this  free  Jet  to  be  highly  unstable  at  a  considerable  distance 
downstream.  Characteristics  associated  with  this  instability 
are  increased  decay  and  spreading  rates.  If  it  were  assumed 
that  the  presence  of  the  flat  plate  and  its  supporting 
structure  had  a  net  stabilizing  influence  on  the  jet  in  this 
axial  range,  the  higher  impingement  pressures  observed  in 
the  downstream  region  could  be  at  least  partially  accounted 
for.  There  is,  however,  no  other  evidence  to  support  such 
an  assumption.  For  the  highly  underexpanded  Jet  p-j/p^  = 

3-57,  the  shift  Is  seen  to  disappear  near  the  end  of  the 
supersonic  portion  x/d^  =  20.  In  the  range  10  <  x/djj  <  18, 
however,  the  apparent  reversal  of  the  shift  may  only  be  the 
result  of  limited  data  taken  in  a  region  where  values  are 
still  quite  sensitive  to  axial  location  because  of  the 
oblique  shock  structure.  Both  curves  in  this  region  are 
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based  on  only  three  data  points.  An  interesting  feature  of 
this  jet  is  the  large  recovery  of  Pitot  pressure  in  the 
core  downstream  of  the  normal  shock  disks. 

Photographic  studies.  Photographs  were  made  of  the  impinge¬ 
ment  flow  on  the  flat  plate  model  by  means  of  schlieren  and 
shadow  techniques.  A  selection  of  the  schlieren  photographs, 
taken  with  a  spark  light  source,  is  reproduced  here  in 
Figures  40  and  4l  for  the  moderately  and  highly  underex¬ 
panded  jets,  respectively.  For  these  pictures,  the  knife 
edge  is  vertical  so  that  density  gradients  in  the  axial 
direction  predominate.  Figure  42  shows  a  series  of  continu¬ 
ous  light  schlieren  pictures  taken  under  similar  conditions 
for  the  highly  underexpanded  jet.  In  the  impingement  region, 
it  is  apparent  that  the  character  of  the  local  structure  as 
revealed  by  the  density  gradients  varies  considerably  with 
impingement  distance.  No  distinct  normal  shock  is  observed 
in  this  region,  although  this  may  be  because  it  is  masked 
by  disturbances  in  the  surrounding  flow  which  are  of  compara¬ 
ble  strength.  As  the  impingement  distance  is  decreased, 
there  is  a  definite  distortion  of  the  jet  structure  upstream. 
For  the  P^/P*,  =  3*57  case  the  axial  location  of  the  Jet 
normal  shock  disk  is  seen  to  start  moving  upstream  for  an 
impingement  distance  between  x/d^  =  6.24  and  5*32.  The 
continuous  light  pictures  show  the  distortion  as  a  sort  of 
"telescoping"  effect  in  the  core.  Possible  changes  in  the 
stability  characteristics  of  these  Jets  due  to  impingement 
cannot  be  deduced  from  these  pictures. 

The  radiation  of  strong  sound  waves  from  the  impinge¬ 
ment  region  is  quite  clear  for  both  jet  strengths.  At 
small  impingement  distances,  there  is  considerable  inter¬ 
action  between  these  waves  and  the  jet  itself . 

Another  series  of  schlieren  pictures  was  taken  in  hopes 
of  revealing  a  visible  correlation  with  the  measured  changes 
in  sign  of  the  radial  pressure  gradient  for  some  of  the 
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cases  cited  earlier.  In  order  to  distinguish  gradients  in 
the  radial  direction,  the  knife  edge  was  held  horizontal. 

With  the  flat  plate  at  an  axial  distance  x/djg  =  1.96, 
the  Jet  pressure  ratio  was  varied  in  small  increments  in  the 
range  1.00  <  Pp/p^  <  5.  This  series  of  pictures  is  shown 
in  Figure  43.  It  is  seen  that  definite  gradients  of  density 
appear  in  a  concentric  pattern  about  the  stagnation  region. 
These  gradients  are  visible  in  every  picture  although  they 
are  relatively  weak  for  the  pressure  ratios  used  in  the 
bulk  of  this  study.  Nevertheless,  they  do  occur  at  radial 
positions  corresponding  to  those  of  the  measured  changes  in 
pressure  gradient.  These  pictures  also  show  the  stand-off 
shock  in  the  impingement  region  very  clearly  for  the  lower 
pressure  ratios.  The  decrease  in  stand-off  distance  with 
increasing  pressure  ratio  is  evident.  Also  of  interest  is 
the  change  observed  when  the  pressure  ratio  is  increased 
from  1.88  to  1.95.  It  is  in  this  range  that  the  normal 
shock  disk  appears  in  the  free  Jet.  In  the  impingement 
case  shown,  the  shock  apparently  forms  at  about  the  same 
pressure  ratio  even  though  its  position  is  shifted  due  to 
the  short  impingement  distance. 

3.2.3.  Evaluation  of  heat  transfer  parameter  (due/dr)r_Q. 

The  stagnation  point  radial  velocity  gradient  (du  /dr)_  n  was 
evaluated  from  the  detailed  stagnation  region  pressure  distri¬ 
butions  given  in  3.2.2.  The  calculations  were  made  by  fit¬ 
ting  the  data  to  Equation  6  at  the  stagnation  point.  In 
order  to  assure  the  best  possible  fit  exactly  at  r  =  0,  the 
pressure  data  were  first  plotted  as  a  function  of  (r/r„)2. 

W 

The  slope  of  the  resulting  curve  was  determined  graphically 
at  r  =  0,  and  the  value  of  (due/dr)r_0  was  then  computed. 
In  most  cases  the  data  fell  close  to  a  straight  line  for  a 
reasonable  distance  from  r  =  0,  so  that  the  slope  evalua¬ 
tion  was  not  difficult.  The  values  computed  in  this  manner 
are  shown  plotted  as  a  function  of  impingement  distance  for 
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each  jet  strength  and  model  shape  In  Figures  44,  45*  46,  and 
47.  Referring  to  the  flat  plate  results,  which  are  more 
detailed,  there  is  a  tendency  for  ( dUg/dr )r=0  to  increase 
at  first  until  a  maximum  is  reached  at  a  point  near  the  end 
of  the  free  jet  core.  Farther  downstream,  a  characteristic 
decay  is  observed.  This  behavior  appears  to  be  a  logical 
consequence  of  the  axial  changes  in  Jet  total  pressure  and 
velocity  profiles,  since  the  radial  gradients  of  these 
quantities  are  smaller  near  the  centerline  in  both  the  core 
and  fully  developed  regions.  In  the  core  of  the  moderately 
underexpanded  jet,  it  is  seen  that  the  values  fluctuate 
sharply  in  a  manner  similar  to  that  observed  for  the  Pitot 
pressure  (Figure  39).  In  the  case  of  the  highly  underexpanded 
jet,  the  negative  value  is  a  result  of  the  reversed  radial 
pressure  gradient  found  to  exist  for  some  distance  downstream 
of  the  normal  shock  disk.  The  existence  of  such  a  gradient 
strongly  suggests  that  there  is  a  separated  region  of  re¬ 
versed  flow  in  the  immediate  vicinity  of  the  usual  stagna¬ 
tion  point.  A  flow  visualization  study  that  seems  to  confirm 

this  condition  is  discussed  in  3*2.4. 

Figures  45,  46,  and  47  show  that  while  the  general  axial 
variation  of  (dug/dr )r=0  is  not  highly  dependent  on  surface 
shape,  the  magnitude  is.  This  dependence  is  made  clear  in 
Figures  48,  49,  and  50,  which  show  the  values  for  all  four 
shapes  plotted  together  for  each  jet  pressure  ratio.  A 
curve  is  drawn  only  for  the  more  detailed  flat  plate  data. 

It  is  seen  that  the  values  for  the  convex  hemisphere  are 
always  higher  and  those  for  the  concave  snapes  always  lower 
than  those  for  the  flat  plate.  This  is,  of  course,  the  same 
as  the  effect  noted  in  the  discussion  of  the  pressure  distri¬ 
butions  themselves. 

Two  methods  of  presenting  measured  values  of  ( due/dr )r=Q 
in  terms  of  measured  jet  behavior  are  given.  The  first 
method,  which  should  be  useful  when  little  or  no  information 
is  available  about  the  decay  and  spread  of  the  impinging  jet. 
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tagnation  point  radial  velocity  gradient  with  impingement  distance 


Figure  48.  Comparison  of  stagnation  point  radial  velocity  gradient  for  a  jet  impinging 

on  several  surface  shapes. 
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Figure  50.  Comparison  of  stagnation  point  radial  velocity  gradient  for  a  jet  impinging 

on  several  surface  shapes . 
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makes  the  stagnation  point  velocity  gradient  nondimensional 
through  the  use  of  d^  and  V^,  i.e. 

\ 

Vxv  dr  /  r=0 

The  second  method,  which  is  more  fundamental  but  which 
requires  a  knowledge  of  the  decay  and  spread  of  the  free 
jet,  makes  the  stagnation  point  velocity  gradient  nondimen¬ 
sional  through  the  use  of  measured  local  free  jet  conditions 
in  the  plane  of  impingement,  i.e. 


The  results  of  the  first  of  these  methods  are  shown  in 
Figures  51-54,  and  those  of  the  second  method  are  shown  in 
Figures  55-58. 

Referring  again  to  the  flat  plate  data,  the  behavior  of 
the  nondimensional  velocity  gradient  based  on  nozzle  exit 
values  is  similar  for  the  two  weakest  jets  (cf.  Figure  44). 
This  similarity  is  indicative  of  the  roughly  equivalent  core 
lengths  and  rates  of  decay  as  well  as  rather  mild  over-all 
effects  due  to  shock  patterns.  For  the  highly  underexpanded 
case,  however,  the  behavior  is  quite  different,  with  the 
effects  of  the  much  longer  core  region  and  strong  shock 
structure  clearly  shown.  Similar  results  are  obtained  for 
the  other  shapes,  but  the  magnitudes  are  different,  as 
expected. 

The  nondimensional  gradient  based  on  local  free  jet 
properties  appears  to  achieve  a  rough  correlation  of  the 
measured  results,  although  the  core  effects  of  the  strongest 
Jet  are  still  in  evidence.  In  general,  the  two  weakest  jets 


Figure  51.  Variation  of  stagnation  point  nondimensional  radial  velocity  gradient  with  impingement 

distance  for  several  pressure  ratios. 


Figure  52.  Variation  of  stagnation  point  nondimens ional  radial;  velocity  gradient  with  impingement 

distance  for  several  jet  pressure  {ratios. 


Figure-  53*  Variation  of  stagnation  point  nondimens ional  radial  velocity  gradient  with  impingement 

distance  for  several  jet  pressure  ratios. 


Figure  54.  Variation  of  stagnation  point  nondimensional  radial  velocity  gradient  with  impingement 

distance  for  several  jet  pressure  ratios. 
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Figure  56.  Variation  of  stagnation  point  nondimens ional  radial  velocity  gradient  v/ith  impingement 

distance  for  several  jet  pressure  ratios. 
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are  well  correlated  about  a  value  (r  ^/Vc)(due/dr )r_Q  ~  1. 1 
for  a  considerable  distance.  Far  downstream,  the  strongest 
jet  also  approaches  this  value.  Aside  from  factors  limiting 
experimental  accuracy,  certain  deviations  may  very  likely 
stem  from  changes  in  the  stability  characteristics  of  the  jet 
due  to  the  presence  of  the  plate.  This  is  particularly  true 
of  the  Pj/p^  =  1.42  jet.  Of  the  other  shapes,  the  best 
correlation  is  shown  for  the  concave  hemisphere.  In  general, 
however,  the  shape  effect  on  magnitude  is  the  same  as  that 
already  noted,  i.e.,  higher  values  for  the  convex  shape  and 
lower  values  for  the  concave  shapes. 

3.2.4.  Visualization  studies  of  stagnation  region  flow.  The 
pressure  distributions  measured  for  impingement  distances 
downstream  of  the  normal  shock  disk  in  the  highly  underex¬ 
panded  jet,  show  peaks  of  maximum  pressure  to  either  side  of 
the  usual  centerline  stagnation  point.  A  typical  example  of 
such  a  distribution  along  with  the  Pitot  and  static  pressure 
profiles  in  the  free  jet  at  the  same  axial  distance  is  shown 
in  Figure  59.  The  pressures  are  all  shown  as  a  fraction  of 
the  free  jet  centerline  Pitot  pressure  so  that  relative 
changes  in  pressure  level  due  to  impingement  may  be  observed. 
The  shape  of  the  free  jet  profile  at  this  point  is,  of  course, 
a  result  of  the  total  pressure  drop  through  the  normal  shock 
disks  upstream.  Even  though  there  has  been  a  substantial 
recovery  on  the  centerline  (see  Figure  39)#  the  outer  region 
of  the  core  retains  its  higher  level  from  the  beginning 
because  of  its  structure  of  relatively  weak  oblique  shocks. 

The  corresponding  impingement  distribution  shows  the  outer 
peak  to  be  displaced  outward,  while  the  central  region  is 
relatively  flat.  At  the  same  time,  the  magnitude  of  the 
outer  peak  is  seen  to  be  even  less  than  that  at  the  center- 
line  of  the  free  jet.  This  condition  suggests  the  existence 
of  a  separated  region  caused  by  the  reversed  pressure  gradient 
to  the  Inside  of  the  outer  peak.  The  flow  pattern  based  on 


Figure  59.  Comparison  of  measured  pressure  distributions  on  flat  plate  and  in 

free  jet. 
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this  idea  is  shown  in  Figure  60.  The  usual  central  stagna¬ 
tion  point  is  transformed  into  a  ring  surrounding  the 
separated  region.  A  pattern  similar  to  this  is  known  to 
exist  under  certain  ground  effects  machines  [27].  A  conse¬ 
quence  of  such  a  flow  in  terms  of  heat  transfer,  of  course, 
would  be  the  existence  of  a  cool  central  spot  surrounded  by 
a  ring  at  stagnation  point  heating  levels. 

Although  confirmation  of  the  phenomenon  by  means  of 
direct  measurement  of  velocity  or  pressure  gradients  off  the 
surface  was  felt  to  be  impossible,  a  method  was  devised  for 
observing  surface  streamlines.  This  method  consisted  of 
impinging  the  jet  for  a  short  time  on  a  layer  of  highly 
viscous  water  pump  grease  applied  to  the  surface  of  the  flat 
plate.  In  order  to  achieve  maximum  visual  contrast  in  the 
resulting  pattern,  the  grease  was  mixed  with  lampblack,  and 
the  plate  was  given  a  smooth  finish  of  white  lacquer.  Since 
it  was  found  that  the  initial  distribution  or  thickness  of  the 
grease  did  not  affect  the  resulting  pattern,  the  application 
method  did  not  present  a  problem.  With  the  grease  applied 
in  a  small  blob,  the  plate  was  covered  with  a  baffle  board 
in  order  to  shield  the  grease  from  transient  phenomena  during 
the  starting  of  the  Jet.  When  the  jet  reached  the  desired 
running  condition,  the  baffle  was  withdrawn  for  a  short  time 
interval  of  about  3  to  5  seconds.  After  replacing  the  baffle, 
the  jet  was  shut  down.  Increasing  the  time  of  exposure  to 
the  jet  was  found  only  to  remove  more  grease  from  the  impinge¬ 
ment  region  and  not  to  change  the  pattern.  Thus,  it  was 
determined  not  only  that  the  phenomenon  was  basically  a 
steady  one,  but  that  it  was  also  quite  repeatable.  A  photo¬ 
graph  of  this  pattern  for  the  highly  underexpanded  Jet 


^Because  of  the  substantial  jet  pressure  force,  the 
baffle  was  made  with  two  raised  edges  that  rested  on  the 
plate  surface  away  from  the  center.  Exposure  of  the  grease 
was  accomplished  by  quickly  sliding  the  baffle  out  of  the  way 
laterally. 


Figure  60, 


Plow  pattern  of  separated  region  (distance  normal 
to  surface  exaggerated). 


Node 


Saddle  point 


ai  surface  streamline  pattern  based  on 
grease  streak  pictures. 
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(p  /p  »  3.57)  at  x/dN  -  5-32  is  shown  in  Figure  62.  The 
behavior  of  surface  streamlines  suggested  by  the  observed 
patterns  is  shown  schematically  in  Figure  6l.  It  is  seen 
that  a  dividing  ring  exists  which  apparently  separates 
regions  of  inward  and  outward  flow.  To  the  inside,  there  is 
a  central  region  of  thick  grease,  while  to  the  outside  an 


alternating  concentric  pattern  of  more  and  less  grease  is 
observed.  Of  particular  interest  is  the  manner  in  which  the 
streaks  curve  near  the  dividing  ring  resulting  in  segmenta¬ 
tion  of  the  central  pattern.  Along  the  dividing  ring,  each 
segment  can  be  identified  by  two  types  of  points  having  the 

topological  character  of  nodes  and  saddle 

points  ^  «  The  number  of  segments  in  the  pattern 

apparently  decreases  with  increasing  impingement  distance. 

For  this  jet,  the  minimum  number  of  three  segments  is 
observed  for  all  distances  beyond  about  x/dN  =4.5  (see 
Figure  63).  Another  factor  which  was  found  to  affect  not 
onlv  the  number  of  segments  at  a  given  distance  but  also 


the  orientation  of  the  pattern,  was  the  jet  nozzle.  Rota¬ 
tion  of  the  nozzle  about  its  axis  between  runs  resulted  in 
an  identical  rotation  of  the  observed  grease  pattern.  It 
was  found  that  tiny  nicks  in  the  nozzle  lip  could  be  corre¬ 
lated  with  at  least  some  of  the  pattern  segments.  Removal 
of  the  larger  nicks  resulted,  for  a  fixed  impingement  dis¬ 
tance,  in  fewer  segments.  The  patterns  shown  in  the  present 
case,  however,  are  those  resulting  from  tests  after  the 
nicks  had  been  removed  from  th;  nozzle.  Under  10-power 
magnification,  no  nicks  were  observed  that  could  be  corre¬ 
lated  with  the  resulting  three  equal  pattern  segments. 

A  comparison  of  schlieren  photograph,  grease  pattern, 
and  surface  pressure  distribution  for  a  typical  case  is 
shown  in  Figure  64.  It  can  be  seen  that  local  pressure  peaks 
are  definitely  related  to  the  observed  patterns,  with  less 
grease  remaining  where  the  pressure  is  higher.  These  local 
pressure  maxima  also  correspond  to  the  dark  regions  of 


Figure  62.  Typical  grease  streak  photograph  of  flow  pattern  due  to 
impingement  on  flat  plate,  p-j/p^  =  3-57;  x/dN  =  5.32. 


I 
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Figure  63.  Grease  streak  patterns  for  impinging  jet;  P^/p^  =  3-57. 
The  scale  divisions  are  tenths  of  an  inch. 


Figure  64.  Typical  jet  impingement  characteristics  for  a  case  in  which  a  normal  shock 

exists  in  the  jet  core.  p../p  =  3.57;  x/d.,  =  5.32. 
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increasing  density  near  the  surface  as  shown  in  the  schlieren 
picture.  In  addition,  the  position  of  the  streamline  divid¬ 
ing  ring,  which  can  be  thought  of  as  a  line  of  zero  radial 
shear,  is  found  not  to  correspond  to  the  local  pressure 
maximum. 


3*2.5.  Momentum  balances  and  interference  effects.  In  order 
to  help  clarify  certain  aspects  of  jet  impingement  flows, 
especially  those  which  may  involve  interference  in  the  en¬ 
trained  flow  field,  it  is  useful  to  consider  qualitatively 
the  balance  of  jet  momentum  fluxes  with  the  pressure  forces 
on  the  impingement  surface.  In  Figure  65,  four  basic  condi¬ 
tions  are  illustrated.  Case  1  shows  a  free  jet  flow  issuing 
from  a  nozzle  of  the  type  used  in  the  present  experiments, 
while  Case  2  shows  this  jet  impinging  on  a  circular  flat 
plate.  Cases  3  and  4  represent  corresponding  flows  but  with 
a  flat  baffle  plate  inserted  at  the  exit  plane  (AB).  For 
convenience,  a  cylindrical  control  volume  is  defined  by  the 
impingement  plate  (surface  CD),  the  baffle  plate  (surface  AB), 
and  the  cylindrical  surface  S.  The  pressure  is  assumed  to 
be  unchanged  across  the  free  boundaries  except  within  the 
confines  of  the  nozzle.  It  is  clear  that  the  total  axial 
momentum  flux  leaving  the  volume  across  surface  CD  in 
Case  1  is  made  up  not  only  of  that  entering  across  AB 
through  the  nozzle  (including  the  effect  of  any  pressure 
difference  across  the  nozzle  exit  plane)  but  also  the  sum  of 
all  the  axial  components  of  entrainment  flow  entering  across 
the  entire  plane  AB  as  well  as  the  surface  S.  Moreover, 
a  complete  force  balance  for  this  system  would  also  have  to 
include  the  proper  components  of  viscous  forces  at  all  the 
boundaries.  For  Case  2,  neglecting  viscous  forces,  the 
balance  is  primarily  between  the  axial  momentum  flux  of  the 
jet  (across  AB  and  S)  and  the  pressure  force  on  the  plate. 
It  should  be  noted,  however,  that  part  of  the  axial  flux 
across  S  is  contained  in  the  wall  jet,  which  because  of  its 
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spread  makes  a  contribution  in  the  negative  axial  direction. 
The  possible  consequences  of  blockage  due  to  a  baffle  plate 
are  clear  in  Cases  3  and  4.  It  is  seen  that  a  large  portion 
of  the  axial  momentum  flux  is  lost  because  of  the  solid 
boundary  at  AB.  The  entrainment  flow  is  thus  primarily 
radial.  Because  this  flow  must  turn  as  it  approaches  the 
jet,  however,  a  pressure  gradient  will  be  established  on 
the  surface  AB  which  can  result  in  an  additional  axial  pres¬ 
sure  force  contribution.  Basically,  Case  4  represents  a 
combination  of  Cases  2  and  3»  Using  the  measurements  of 
free  jet  velocity  profiles  and  over-all  flat  plate  pressure 
distributions,  a  comparison  between  jet  axial  momentum  flux 
and  plate  pressure  force  has  been  made  for  the  Pt/p^  =  1*42 
jet.  Each  of  these  quantities  is  plotted  as  a  function  of 
axial  distance  in  Figure  66.  In  addition,  the  free  jet  mass 
flux,  the  radial  boundary  of  the  nozzle  mass  flux  rj5[N/rN 
(i.e.,  the  dividing  streamline  between  nozzle  and  entrained 
flow),  and  the  over-all  jet  spreading  parameter  r^/r^  are 
given  in  Figure  67.  The  behavior  of  the  jet  momentum  flux 
and  plate  force  can  be  explained  in  terms  of  that  portion  of 
the  over-all  axial  flux  actually  measured  at  each  station. 

Near  the  nozzle,  of  course,  the  plate  force  represents 
practically  all  of  the  axial  flux  due  to  the  nozzle  as  well 
as  entrainment.  It  is  seen  that  the  jet  alone  accounts  for 
only  about  one -half  of  this  total  force.  The  initial  decrease 
of  each  parameter  with  axial  distance  is  believed  to  be 
primarily  due  to  pressure  gradients  which  are  known  to  exist 
across  the  boundary  AB  near  the  nozzle,  but  which  were 
not  evaluated.  Farther  downstream, the  portion  of  the  total 
control  volume  cross-section  included  in  the  free  jet  measure¬ 
ments  increases  because  of  jet  spreading  so  that  a  larger 
portion  of  the  total  axial  flux  is  accounted  for.  At  the 
same  time,  since  the  jet  diameter  has  grown  to  be  nearly 
that  of  the  plate,  some  of  the  jet  flow  is  deflected  around 
the  plate  without  contributing  to  the  force.  This  drop-off 


> 
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Figure  67.  Measured  total  mass  flux  in  free  jet  and  spread 
of  nozzle  and  total  mass  flux. 
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is  slight  within  the  range  of  these  measurements. 

The  measured  jet  mass  flux,  as  shown  in  Figure  67,  is 
seen  to  increase  more  rapidly  in  the  downstream  portion.  The 
possible  existence  of  instabilities  in  this  region  has  already 
been  discussed.  Because  of  probable  measurement  uncertain¬ 
ties  in  such  a  region,  therefore,  it  is  not  clear  how  much 
of  the  increased  flux  is  actually  due  to  increased  entrain¬ 
ment.  Nevertheless,  the  total  mass  flux  at  x/dN  ss  40  is 
seen  to  be  approximately  l4  times  that  of  the  Jet  alone.  In 
order  to  illustrate  the  interference  effect  of  a  baffle  plate 
at  surface  AB,  several  impingement  conditions  were  rerun 
both  with  and  without  such  a  baffle  plate.  A  comparison  of 
impingement  pressure  distributions  resulting  from  this  test 
is  shown  in  Figure  68.  The  change  in  distribution  in  this 
particular  case  is  characterized  by  a  decreased  pressure 
near  the  center  and  an  increase  in  the  outer  portions. 

While  some  of  the  few  cases  tested  exhibited  similar  changes, 
others  did  not.  On  the  whole,  the  relationship  between  jet 
and  impingement  behavior  and  interference  effects  is  extreme¬ 
ly  complex  and  the  present  simple  tests  are  intended  to 
provide  only  some  idea  of  the  changes  that  can  result  from 
changes  in  geometry. 

3«3.  Discussion. 

The  experimental  determination  of  stagnation  point 
heat  transfer  parameters  for  normal  jet  impingement  has  been 
described  in  detail.  It  has  been  found  that  by  making  use 
of  local  free  jet  characteristics  a  reasonable  correlation 
of  the  stagnation  point  velocity  gradient  (due/dr)r=Q  is 
achieved  for  a  number  of  cases  of  Jet  strength  and  impinge¬ 
ment  distance.  For  impingement  on  a  flat  plate  it  was  found 
that  ( r . 5/^ ) ( due/dr ) r _q  fell  within  20  per  cent  of  a 
value  of  1,1  at  axial  distances  between  10  and  40  nozzle 
diameters  downstream  for  a  subsonic  (p^/p^  =  .800)  Jet  and 
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a  moderately  underexpanded  jet  (p-j/p^  =  1.42).  For  a  highly 
underexpanded  jet  (p-j/p^  =3.57)  a  value  of  1.1  was 
approached  within  20  per  cent  downstream  of  30  nozzle 
diameters.  For  a  convex  hemisphere,  generally  higher  corre¬ 
lation  values  were  found.  These  values  showed  an  upward  trend 
with  increased  axial  distance.  Such  a  trend  might  be  expected 
since  far  downstream  the  axial  drop-off  of  the  gradient 
(due/dr)r_Q  itself  is  less  rapid  than  it  is  for  the  flat 
plate.  While  the  concave  shapes  exhibit  lower  values  than 
for  the  flat  plate,  data  do  not  clearly  indicate  the  downward 
trend  which  might  be  expected  on  similar  grounds.  The 
correlation  values  can  be  related  to  the  geometric  conditions 
in  terms  of  the  free  jet  half-velocity  radius  r  and  the 
model  surface  radius  of  curvature  rg.  With  ,:r^  defined 
negatively  for  the  concave  hemisphere,  the  values  for  the 
fully  developed  region  of  each  jet  are  shown  in  Figure  69. 

In  the  range  -0.6  <  r  <  0.6  ,  the  surface  curvature 

effect  can  be  approximated  by  the  linear  relation 


1.13  +  1.08 


In  addition  to  these  basic  measurements  of  (due/dr)r_0  , 
several  other  measurements  and  photographic  studies  have 
been  made  which  help  to  illustrate  certain  features  of 
impingement  flows,  particularly  of  highly  underexpanded  jets. 
It  has  been  found  that  under  conditions  of  impingement  out 
to  about  8  nozzle  diameters  for  a  jet  with  p-^/p^  =  3.57,  a 
separated  region  may  exist  in  the  vicinity  of  the  stagnation 
point  witn  the  result  that  maximum  heat  transfer  may  occur  in 
a  ring  surrounding  the  central  region.  Some  unusual  features 
of  this  flow  have  been  investigated  by  comparing  impingement 
pressure  distributions,  surface  streamline  patterns,  and 
schlieren  photographs.  Of  special  interest  has  been  the 
manner  in  which  the  separated  flow  region  becomes  segmented 
along  the  ring  dividing  the  inward  from  the  outward  flow  on 
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the  surface.  It  has  also  been  found  that  the  maximum 
impingement  pressure  occurs  outside  of  this  ring. 

The  problem  of  interference  due  to  obstructions  in 
the  flow  field  surrounding  the  jet  proper  has  been  discussed 
briefly  in  a  qualitative  way.  In  particular,  it  has  been 
stressed  that  the  presence  of  such  obstructions  can  result 
in  substantial  changes  in  the  impingement  pressure  distri¬ 
bution.  An  example  of  this  effect  has  been  given  in  which 
interference  was  introduced  by  placing  a  baffle  plate  in 
the  plane  of  the  nozzle  exit. 
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4.  CONCLUSIONS 
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This  study  was  performed  in  order  to  obtain  a  general 
knowledge  of  the  behavior  of  free  jets  when  they  impinge 
normally  on  a  flat  surface.  In  particular,  the  general 
character  of  the  stagnation  point  heat  transfer  parameter 
(due/dr)r_Q  and  its  dependence  on  shock  structure  and  other 
features  of  the  free  Jet  were  studied  in  detail. 

The  purpose  of  this  investigation  was  to  provide  a 
general  guide  to  the  study  of  such  matters  at  high  enthalpy 
where  the  effects  of  density  variations  are  of  importance. 

It  was  felt  that  by  knowing  the  general  features  of  free  Jet 
impingement  for  low  temperature  constant  density  flows,  a 
great  deal  of  labor  and  expensive  high  temperature  testing 
could  be  eliminated  and  a  limited  but  essentially  definitive 
test  program  for  the  high  enthalpy  problem  evolved. 

The  results  of  the  low  temperature  studies  presented 
herein  are  thought  to  be  sufficient  to  permit  the  order  of 
magnitude  of  the  stagnation  point  heat  transfer  parameter 
(due/dr)r_0  to  be  estimated  for  the  case  of  high  enthalpy 
Jets  in  a  number  of  important  practical  applications. 

The  results  of  these  low  temperature  studies  of 
(due/dr)r_0  were  presented  in  two  nondimensional  forms.  In 
the  first  form  (due/dr)r_0  was  made  dimensionless  using  the 
Jet  exit  diameter  and  jet  exit  velocity,  i.e. 

V1  l  dr  /r=0 

This  form  is  useful  for  quick  estimates  of  (dup/dr)r_Q  when 
information  is  not  available  concerning  the  decay  of  the  free 
Jet  under  consideration.  The  other  method  of  presenting  the 
data  in  nondimensional  form  attempts  to  correlate  the  stagna¬ 
tion  poini-  heat  transfer  parameter  with  local  conditions  in 
the  free  jet  at  the  plane  of  impingement,  i.e. 
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The  data  show  that  such  a  correlation  is  useful  in  the  region 
of  fully  developed  free  jet  flow.  In  particular,  for  the 
case  of  normal  impingement  on  a  flat  plate 


The  effect  of  curvature  of  the  impingement  surface  on 
the  stagnation  point  heat  transfer  parameter  was  also  investi¬ 
gated.  Typically,  for  fully  developed  jets  when  the  ratio  of 
the  half  width  of  the  free  jet  r  _  at  the  plane  of 
impingement  to  the  radius  of  curvature  r^  of  the  impingement 
surface  falls  in  the  range  -0.6  <r^/rs<  0.6  (the  plus  sign 
refers  to  a  convex  surface  and  the  minus  sign  to  a  concave 
surface),  the  effect  of  surface  curvature  on  (r  ^Vc )  ( due/dr )  r_Q 
is  given  approximately  by  the  expression 


1.13  +  1.08 

rS 


In  order  to  use  this  second  method  of  nondimensionalizing 
the  stagnation  point  heat  transfer  parameter,  it  is  necessary 
to  compute  the  properties  of  a  given  jet  downstream  of  its 
point  of  issuance.  This  problem  was  also  investigated. 
Comparisons  were  made  between  measured  and  computed  values 
of  r  c  and  V  at  various  positions  along  the  jet.  Computed 
values  were  found  by  applying  Warren’s  momentum  integral  method 
to  all  the  cases  studied  including  those  for  which  the  jet 
was  underexpanded.  It  was  found  that  the  predicted  decay 
rates  based  on  this  method  agreed  with  the  measured  values 
quite  well  in  the  subsonic  range.  For  moderately  underexpanded 
jets,  however,  the  agreement  was  no  better  than  about  60  per 
cent.  This  is  believed  to  be  due  primarily  to  the  effects  of 
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a  large  scale  Instability  or  "flapping"  observed  for  these 
jets.  The  agreement  was  much  better  in  the  highly  under¬ 
expanded  range,  although  a  steady  drop-off  in  agreement  with 
increasing  jet  strength  and  axial  location  was  noted.  Since 
the  analytical  method  was  based  on  a  turbulent  exchange 
coefficient  that  was  assumed  to  be  a  function  only  of  Initial 
jet  Mach  number,  it  is  felt  that  such  a  drop-off  might  be 
expected  as  the  jet  decays  to  subsonic  velocities.  These 
data  would  seem  to  indicate,  therefore,  that  in  the  absence  of 
jet  "flapping",  the  core  shock  structure  does  not  have  a 
first-order  effect  in  determining  the  decay  rate  of  free  jets. 
It  appears  that  the  disagreement  is  due  primarily  to  the 
limitation  imposed  by  assuming  a  constant  exchange  coefficient 
independent  of  position  along  the  jet.  Only  in  cases  where 
shocks  contribute  to  jet  instability  does  it  appear  that  they 
may  influence  over-all  decay  processes. 
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APPENDIX  I 


Summary  of  free  jet  cases  for  which 
velocity  profiles  were  determined. 
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In  addition  to  the  above  listed  cases,  axial  decay 
measurements  alone  were  made  at  several  of  the  same  locations 
for  the  following  jet  pressure  ratios: 

Po/Psc  =  *244j  pl/p~  =  2’16 

.183  2.88 

.121  4.36 

.104  5.08 
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Pressure  measurement  data  and  computed  velocity  profiles 
for  all  of  the  free  jet  cases  listed  in  Appendix  I  are 
included  in  the  following  pages. 

Summary  plots  of  the  radial  spread  and  axial  decay 
characteristics  of  these  Jets  are  also  included* 


Notes  For  ease  in  making  comparisons,  all  measured 
pressure  data  have  been  converted  to  the  same 
units --inches  of  alcohol  (s.g.  =  .820)  gauge. 
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Figure  AII-10 


Figure  AII-14.  Radial  Spread  of  Underexpanded  Free  Jet 


Figure  All -15.  Axial  Decay  of  Centerline  Velocity  in  Subsonic  Free  Jet 


appendix  III 


Summary  of  impingement  cases  for  which  radial  velocity 
gradient  was  determined  from  detailed  stagnation  point 
pressure  distribution.  This  matrix  was  repeated  for 
each  of  the  four  impingement  models. 
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Additional  cases  were  as  follows: 


For  the 

flat  plate  model. 

P«/psc 

=  .961, 

*/aN 

=  7.32 

P«/Psc 

=  .462, 
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=  7.32 
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Measurements  oi  over-all  surface  pressure  distributions 
were  made  for  all  four  models  for  the  cases  listed  in 
the  main  matrix  above. 


3.  Stagnation  region  pressure  distributions  were  also 
measured  under  the  following  conditions  with  the  flat 
plate  model  (these  data  are  not  all  treated  specifically 
in  this  report): 


p«/psc  varied  ln  small  increments,  x/dM  =  1.96, 
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23.5, 
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Po/Psc  *  -3TS, 
P«/P°c  - 


x/d^  varied  in  small  increments 
x/d^j  varied  in  small  increments 


APPENDIX  IV 


Pressure  distribution  data  for  each  basic  test  condi¬ 
tion.  Over-all  surface  pressures  are  plotted  as  a  percentage 
of  jet  settling  chamber  stagnation  pressure. 
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Figure  AIV-4.  Impingement  pressure 
distributions  for  several 
impingement  distances . 
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Figure  AIV-6.  Impingement  pressure 
distributions  for  several 
impingement  distances. 
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Figure  AIV-7.  Impingement  pressure 
distributions  for  several 
impingement  distances. 
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Figure  AIV-8.  Impingement  pressure 
distributions  for  several 
impingement  distances . 
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Figure  AIV-11.  Impingement  pressure 
distributions  for  several 
impingement  distances. 
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